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Summary 


A multivariable control design method based on constrained parameter optimization 
has been applied to the design of a multiloop aircraft flight control system. Specifi- 
cally, the design method is applied to 

• Direct synthesis of a multivariable “inner-loop” feedback control system based 
on total energy control principles, 

• Synthesis of speed/altitude-hold designs as “ outer-loop” feedback/feedforvvard 
control systems around the above inner loop, 

• Direct synthesis of a combined “inner-loop”and “outer- loop” multivariable con- 
trol system. 

The design procedure offers a direct and structured approach for the determination of 
a set of controller gains that meet design specifications in closed-loop stability, com- 
mand tracking performance, disturbance rejection, and limits on control activities. 
The presented approach may be applied to a broader class of multiloop flight control 
systems. Direct tradeoffs between many real design goals are rendered systematic 
by this method following careful problem formulation of the design objectives and 
constraints. Performance characteristics of the optimized design have been improved 
over the current autopilot design on the B737-100 Transport Systems Research Vehi- 
cle (TSRV) at the landing approach and cruise flight conditions; particularly in the 
areas of closed-loop damping, command responses and control activity in the presence 
of turbulence. 
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Chapter 1 

INTRODUCTION 


Recent development in multivariable control design techniques have been focused 
on the improvement of controller performance and robustness to uncertainties in the 
plant model. Modern control system design methods such as LQG/LTR, H°°, and 
^i-synthesis can provide controllers with high levels of performance and robustness 
[Refs. 1-3]. However, the controllers obtained with these techniques are usually of 
high order. The number of controller states is usually greater than or equal to the 
number of states in the plant model. Controller order must subsequently be reduced 
for the final design implementation. In addition, these design methods do not pro- 
vide direct tradeoffs between many real design performance and robustness measures. 
Alternative methods based on output feedback are available [Refs.4-6] to synthesize 
controllers of low dimensionality. In particular, the design algorithm developed in 
Reference 6 provides a convenient framework for designing optimal low-order con- 
trollers. Recent extension of this method [Ref.7] and the use of constrained nonlinear 
optimization [Ref.8] enable designers to address conveniently those requirements, for 
example, eigenvalue and covariance constraints, that are not easily expressible in 
terms of an objective function based on quadratic performance indices. 

Traditionally, classical design procedures based on single-loop closure and root- 
locus have been applied to the design of multiloop control systems. Designs obtained 
with conventional methods can have minimally achieved performance and inadequate 
robustness to plant model uncertainties. Tradeoffs between performance and robust- 
ness requirements are far from simple and can be time-consuming; especially when 
the designer lacks appropriate past design experiences or when he is faced with a non- 
conventional design case. Furthermore, the classical design procedures often overlook, 
due to their single-input/single-output nature, the intrinsic multivariable aspects of 
the design problem. As a result, the controllers may not contain appropriate cross- 
feed among individual feedback loops, and thereby limiting the maximally achievable 
performance. 
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Complexities in present automatic flight control systems (AFCS) for commercial 
transport aircraft are generally associated with the lack of suitable control system 
design integration, i.e crossfeed among different control paths. Recent work by 
Lambregts [Refs .9- 1 1 ] on the NASA B737-100 Transport System Research Vehicle 
(TSRV), in improving the operation of the AFCS, has led to the development of a 
total energy control system (TECS) for an integrated autothrottle/autopilot design. 
The design philosophy incorporates fundamental aspects of the aircraft dynamics in 
the formulation of a multiloop controller structure. This innovative design approach 
has led to a controller structure that contains key cross-coupling between the elevator 
and throttle control loops for improved performance. Selection of the feedback gains 
was still performed using classical design procedures [Refs. 9-14]. To satisfy multiple 
design requirements such as closed-loop stability, control and command bandwidths, 
control activities in the presence of turbulence, design robustness to modeling uncer- 
tainties, the one-loop-at-a-time iterative procedure can be overly time-consuming. 

Initial design and evaluation of the TECS concept were performed on the NASA 
B-737-100 Transport Systems Research Vehicle (TSRV) [Refs. 12- 14], Selection of 
the feedback gains was still a labor-intensive process even though the philosophy of 
TECS has inherently reduced the number of design parameters in the control-law 
structure. Again classical SISO design procedures without the benefits of past design 
experiences provides little insight into the design tradeoff; particularly when the num- 
ber of constraints and design variables become large. Resulting “inner-loop” designs 
may possess unnecessarily high bandwidths in the control paths. Furthermore, gain 
scheduling to different flight conditions and robustness to plant model uncertainties 
pose additional complexities to the design problem; hence making the SISO design 
procedure lesser attractive than a direct procedure based on constrained parameter 
optimization. 

Multivariable control design techniques could potentially provide a systematic ap- 
proach to the solution for a set of feedback gains that achieve simultaneously de- 
sign performance and robustness to plant model uncertainties. One such approach 
is provided by the unified design algorithm for robust low-order controllers [Ref.6] 
using constrained parameter optimization as fully implemented in the computer pro- 
gram SANDY. The design method has found numerous applications in flight controls 
[Refs. 15-17] and control of flexible mechanical systems [Ref.18]. The work conducted 
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in this study is to demonstrate the application of the design algorithm to the syn- 
thesis of an integrated autothrottle/autopilot ’’inner” structure based on the TECS 
design philosophy. In addition, we perform autopilot designs for speed and altitude 
holds based on an existing optimized inner-loop and compare the results with those 
obtained from simultaneous redesign of both the inner (i.e TECS) and outer-loops. 

The main objective of this work is the formulation of an optimal control problem 
addressing the same design objectives as those achieved under classical design ap- 
proaches; in this case the design solution is based on direct constrained optimization. 
Within this formulation, a systematic approach is developed allowing the design- 
ers to achieve multiple design objectives such as damping of system modes, shaping 
of loop frequency responses, disturbance rejection, command tracking, and robust- 
ness to modeling uncertainties. Using this design technique, feedback gains have 
been obtained for the integrated autothrottle/autopilot design on the TSRV vehicle 
at two flight conditions. Linear analysis indicates that the design procedure using 
constrained optimization yields comparable and improved results in closed-loop sta- 
bility, aircraft responses and control activities to turbulence, and robustness in terms 
of single-input/single-output gain and phase margins. 

This study is divided into three parts: 

1. Synthesis of an integrated autothrottle/autopilot “inner”-loop structure based 
on the TECS design philosophy, 

2. Synthesis of airspeed-hold and altitude-hold “outer”-loops with the previously 
designed “inner”-loop, 

3. Synthesis of an integrated autothrottle/autopilot using simultaneously design 
parameters from both the “inner” and “outer”-loop structures. 

A systematic approach to multivariable controller design is developed allowing the 
designer to perform tradeoff among different performance objectives such as damping 
of system modes, command tracking, shaping of loop frequency responses, disturbance 
rejection, and robustness to modeling uncertainties. Direct tradeoffs between design 
criteria are conducted in a systematic manner, enabling the designer to explore the 
full potentials of the control system. The TECS control system designed under this 
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approach is found to be sufficiently insensitive to variations in the plant model, while 
disturbance rejection requirements to clear-air turbulence and responses to flight path, 
acceleration, airspeed and altitude commands are simultaneously improved. 



Chapter 2 

INTEGRATED AFCS DESIGN 

2.1 Review of AFCS Designs 

Development of automatic flight control systems (AFCS) for modern transport air- 
craft has resulted from the need to reduce pilot workloads and, at the same time, 
improve performance and fuel efficiency. The first flight control systems consisted of 
a pitch attitude hold system through simple feedback of pitch attitude and pitch rate 
to the elevator control. Speed-hold and altitude-hold modes followed soon afterwards 
as extentions of the basic flight control system. Design of new control modes has often 
proceeded with little or no consideration of existing modes in the AFCS. As a result, 
there is numerous replication of the basic flight-path and speed control functions in 
the navigation and control computer systems [Ref.10]. 

In recent years, the need to resolve fundamental operational and performance 
deficiencies of current AFCS designs has become apparent. Extensive development 
programs [Ref. 13] have shown that the conventional AFCS design has reached a 
fundamental limit such that further improvements cannot be easily obtained within 
the existing system architecture. 

The fundamental limitation resides in the traditional single-loop design approach 
where the throttle and elevator commands are developed separately without regard 
to cross-coupling effects in the longitudinal aircraft dynamics. This type of single- 
loop approach is evident in the underlying structure of conventional AFCS designs, 
Figure 2.1. Flight-path control by the autopilot is achieved through feedback to the 
elevator while the autothrottle independently controls speed through the throttles. 

The basic AFCS structure has serious design deficiencies. In the longitudinal air- 
craft dynamics both the elevator and throttle controls produce responses in flight path 
and speed. This effect may lead to adverse cross-coupling between the autothrottle 
and autopilot when both feedback loops are closed. Performance and stability of the 
autothrottle loop may significantly degrade in the presence of the autopilot loop. In 
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fact, for large flight-path angle commands, the much larger control authority and 
bandwidth of the elevator can easily out strip the throttle’s ability to maintain the 
commanded speed. 

As a result, in 1979 NASA directed the Boeing Company to begin conceptual 
development of an integrated flight path and speed control system incorporating many 
existing design objectives while avoiding duplication of control functions. Numerous 
types of “energy compensation” techniques were investigated and led to the design 
and implementation of the so-called total energy control system (TECS) [Refs. 9- 
11]. The total energy control concept provides an effective means of dealing with 
aircraft operational requirements, control nonlinearities and performance limitations 
[Refs. 12- 14]. 

Key features of TECS are 

• Integration of vertical control modes into a single control concept designed 
around a fixed “inner-loop” structure, 

• Use of cross-feed paths between flight path and acceleration to both the elevator 
and the throttles, hence formulation of a multiloop control design structure, 

• Control strategies based on total energy rate and total energy distribution rate, 

• Command paths for control of flight path and longitudinal acceleration. 

These contrast with the conventional approach to AFCS design involving separate 
pitch autopilot and autothrottle designs; in which the elevator is directed to control 
flight path while airspeed regulation is acquired separately through the throttles. 

2.2 Design Objectives 

Developmental studies have been directed toward integrating numerous functions of 
the flight control system into a single control concept. Several important considera- 
tions in the redesign of the AFCS are 

• Integration of all vertical control modes into a single control-law concept. The 
system should be designed around a fixed inner-loop configuration with outer- 
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loop control modes generating signals compatible with the inner- loop command 
inputs. 

• Minimization of transient errors due to cross coupling between commanded in- 
puts in the inner-loop structure. The control system should produce a “coor- 
dinated” response between the throttle and the elevator controls to command 
inputs. 

• Elimination of functional overlaps in the AFCS design. 

• Protection against exceeding aircraft performance and structural limits by pro- 
viding an intelligent hierarchy in control system modes. 

• Provision of overshoot-free responses to step command inputs. 

• Constraint of closed-loop damping of dominant poles to an acceptable level. 

• Minimization of control activities to turbulence in accordance with meeting 
flight-path and speed tracking requirements at each flight condition. 

2.3 Review of the Total Energy Control Concept 

Numerous variations of “energy compensation” techniques have been investigated 
for improving the simultaneous operation of the autothrottle and autopilot designs 
[Ref. 1 3]. Energy compensation becomes the foundation of an integrated AFCS “inner- 
loop” structure in TECS. One formulation of the energy compensation technique is 
to develop the engine thrust command as a function of the aircraft total energy 
rate, while the elevator command is expressed as a function of the distribution rate 
between kinetic and potential energies. Key equations governing the total energy 
control concept were developed by Lambregts [Refs. 9- 11]. They are included here for 
completeness and to define the structure of the control system used in the constrained 
parameter optimization that follows. 

The total energy E(t) of the aircraft treated as a point mass is given by 

E(t) = -mV(t) 2 + mgh(t) 


(2.1) 
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where 

m — mass of aircraft 

V(t ) = aircraft total velocity along the flight path 

h(t ) = altitude 

g = gravitational acceleration 

The total energy rate E(t) is found by differentiating E(t) in equation (2.1) with 
respect to time as, 

~ rngV{t) + 7(0^ (2-2) 

for small flight-path angle 7 (t) in radians and assuming that the aircraft mass m is 
constant or slowly varying. If we write the equation for the thrust required from the 
equations of motion along the flight path, we have 

T,„(t) = mg fnil + 7 (f)j + D(t) (2.3) 

Assuming that initial thrust is trimmed against drag and the variations in drag with 
time are generally slow, then it follows that the short-term thrust requirement for 
a desired level of total energy rate is obtained from the acceleration and flight-path 
angle quantities as 

Ar„,(() « ||| « mg fM + 7 (()j (2.4) 

This implies that the aircraft total energy can be regulated directly using thrust 
control. However the distribution between kinetic and potential energies cannot be 
effectively controlled with the throttles. To distribute the total energy rate between 
flight path and acceleration as desired, elevator control must be used. The elevator 
primarily provides control of angle-of-attack and thus lift while contributing little to 
the aircraft total drag, thereby having little effect upon the total energy rate. 

Using proportional and integral control on the total energy rate and energy rate 
distribution, we develop command strategies for both the thrust and the elevator 
controls as 

S r ,(s) = mg (k tp + — p) ( 77 ^ + 7.(3)) (2.5) 

= (k ep + ^i) - 7,w) 



( 2 . 6 ) 
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where 7 { (s) = 7(5) — 7 c (-s) and K(s) = V(s) — K(s) represent the errors of flight-path 
angle 7 and acceleration V from the commanded values 7 c and V c respectively. The 
parameters Ktp, Kti, Kep, and Kei are the proportional and integral feedback gains 
to thrust and elevator controls respectively. The basic implementation of this system 
is shown in Figure 2.2. The engines can be made to produce the required thrust by 
converting the total thrust command into an equivalent engine pressure ratio (EPR) 
command and closing an EPR feedback loop around the engine [Ref. 1 1]. Notice that 
proportional paths on the flight path 7 C and acceleration V c commands have been 
left out of the current implementation. It was found that direct feedthrough of the 
commanded inputs to the controls would produce undesireable overshoots. However 
such a problem can be resolved with proper feedforward control as evident in this 
study. 

2.4 Analysis of the Current Total Energy Control Structure 

A careful study of the TECS structure and its characteristics is crucial to the for- 
mulation of an optimal control design procedure. The NASA TSRV B-737 research 
aircraft was used to define linear longitudinal aircraft dynamic models at two flight 
conditions (Appendix E). The chosen flight conditions are listed in Table 2.1. The 
first flight condition FLTl is a typical approach to landing configuration, and the 
second flight condition FLT2 is a cruise condition. 

Table 2.2 shows the inner-loop and pitch-damper gain selection for the current de- 
sign implementation on the NASA TSRV [Refs. 12-14]. The gains Kg and K q provide 
feedback of pitch attitude and pitch rate for stabilization of the aircraft longitudinal 
rigid-body modes. Notice that the value of Ktp was set to zero, thus eliminating 
the proportional feedback path to the throttles. Classical analysis indicates that this 
feedback gain is not neccessary for closed-loop stability. Furthermore, by eliminat- 
ing direct feedback to the throttle, the broken-loop control bandwidth and throttle 
activity to turbulence are reduced. Actually this restriction on Kjp is not needed 
as demonstrated in the optimal TECS design. The additional degree-of- freedom pro- 
vided by the gain Ktp can be used to further improve stability of the phugoid and 
short-period modes, and to achieve better command tracking performance without 
increasing control activity response to turbulence. 
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The integral gains Kji and Kei were chosen to be approximately equal so as to 
obtain equal bandwidths in both the 7 e and V t /g feedback paths. Gain scheduling of 
the control system to different flight conditions is achieved through the gains I\cas 
and Kqw in the elevator and throttle control loops respectively (Figure 2.2). The 
gain schedules are a function of calibrated airspeed (CAS) in fps and gross weight in 
pounds and they are given as follows, 


Kcas 

Kgw 


r 200 \ 2 

\CASJ 

mg 


(2.7) 

( 2 . 8 ) 


With these gain schedules the control-law when implemented will clearly be nonlinear 
and time-varying. However in the neighborhood of a design condition, one can still 
examine the responses in terms of those from a linear time-invariant system. 

The TECS structure as shown in Figure 2.2 introduces an uncontrollable pole at 
the origin in the closed-loop system. This results from the additional integration of 
the aircraft acceleration variable V(t) in the control system. This means that the 
TECS design provides only one integral control action effective on the flight-path 
response and not on velocity. Integral control of the velocity variable can only be 
obtained through feedback of the velocity error to the inner-loop commands 7 C and 
VJg, as found in an outer-loop control of a speed-hold type autopilot. Design of 
outer-loops for different autopilot modes (i.e speed-hold and/or altitude holds) will 
be discussed in chapters 7 and 8. In this section, we are primarily concerned with the 
redesign of the TECS “inner”-loop gains for improved performance and robustness 
using a procedure based on nonlinear constrained optimization. 

Figure 2.3 shows frequency responses from flight-path angle and acceleration com- 
mands to the respective outputs for the closed loop system at flight condition FLTl. 
Note that at zero frequency the transfer function of V(s)/1/ C (s) has a gain less than 
one. Consequently the controller will not maintain the commanded acceleration in 
steady state. The TECS architecture simply cannot provide a constantly increasing 
thrust command (i.e. a type 2 system in the acceleration command path) when the 
error quantities V ( (t) and 7 f (<) become zero, i.e. when V"(t) = V c (t ) and q(t) = 7 C (0- 



Table 2.1: TSRV B-737 Flight Conditions 



Altitude 

(ft) 

Weight 

(lbs) 

CAS 

(kts) 

Flaps 

(deg) 

FPA 

(deg) 

Gear 

FLT1 

1,500 

80,000 

120 

40 

-3 

down 

FLT2 

25,000 

80,000 

450 

0 

0 

up 


Table 2.2: Current TECS Inner-Loop and Pitch Damper Gains 


Kti 

I\ei 

Kjp 

Kep 

I< e 

K q 

0.4 

2.52 

0.0 

3.36 

6.0 

4.0 



Figure 2.1: Conventional AFCS Autothrottle/ Autopilot 
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Figure 2.2: Basic Total Energy Control System 
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Frequency (Rad/Sec) 

Figure 2.3: TECS Command Loop Frequency Responses 



Chapter 3 

OPTIMAL CONTROL DESIGN ALGORITHM USING 
CONSTRAINED PARAMETER OPTIMIZATION 


Development of techniques for multivariable control-law synthesis is aimed at pro- 
ducing controllers with desired performance in command tracking and disturbance 
rejection along with adequate robustness to plant uncertainties. Numerous multi- 
variable control design techniques are available such as standard linear quadratic 
gaussian (LQG), linear quadratic gaussian with loop transfer recovery (LQG/LTR), 
//“-optimization and /^-synthesis [Refs. 1-3]. However, simple design procedures for 
low-order controllers addressing direct tradeoff between performance and robustness 
are still lacking. The procedure based on nonlinear constrained optimization as imple- 
mented in the computer program SANDY [Refs. 6, 7, 15, 21] seems to offer a practical 
alternative amidst the abundance of multivariable robust control design schemes. In 
this section, we describe the formulation of models for control-law synthesis, and the 
definition of various performance indices and constraints used in the design algorithm. 

3.1 Formulation of Models for Control-Law Synthesis 

Aircraft dynamic models for linear control-law synthesis are usually obtained from the 
linearization of a full nonlinear model in real-time simulator about some predefined 
equilibrium points. A linearized plant at the i th condition is modelled as a continuous 
linear time-invariant system of the form 


*;(o 

= A' p x' p (t) + B' p u' p (t) + I>J,(<) 

(3.1) 

yl(t) 

= c;xj,(0 + D‘y p (t) + n>’(0 

(3.2) 

4(0 

= c p x' p (t) + £>x(f) + a >'(() 

(3.3) 


where the matrices A p , B p , IT , C\, D\ , Cl' a , C p , D p and fi], contain the plant dynamics, 
actuator and sensor models, disturbance and command generator filter models, etc., 
(i = l,...,JV p ) and N p is the total number of design conditions. The output vector 
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y',(t) represents the sensor variables of the i th plant condition. The vector y' s (t ) in 
equation (3.2) contains not only the measured feedback quantities such as pitch rate, 
pitch attitude, airspeed, etc..., but also the tracked commands (Appendix B) which 
are outputs of a command generator model. This model is often imbedded as part of 
the augmented plant description. For example, y 3 (t) = {</,#, 7, V, 7 C , V^}. The output 
vector y' p (t) contains those outputs used in the performance index of section 3.2 . 

The design method based on nonlinear constrained optimization requires a-priori 
the setup of a predefined controller structure for both the feedback and feedforward 
control systems. The form of the controller is quite general. The complete arbi- 
trariness allowed in the controller structure means that designers will have the direct 
responsibility to: 

1. Setup a control-law with adequate number degrees of freedom so that the design 
algorithm will converge to a solution that meets the desired performance goals 
and the specified set of design constraints and, 

2. Avoid overspecifying the design parameter set (especially in the case of a nonmin- 
imal state realization of the controller). This could lead to nonunique solutions 
or to the breakdown of the search algorithm. 

A typical feedback/feedforward controller has the form 

x c (<) = A c x c {t) + D c u c {t) (3.4) 

y c (t) = C c x c (t) + D c u c (t) (3.5) 

where the state matrices /l c , B c , C c and D c describe the form of dynamic compen- 
sation (i.e. proportional, integral, derivative, leads, lags, etc.) applied to the input 
vector u c (t). The controller output vector y c (t) contains the feedback controls and 
those output variables that are expressed as linear combination of controller states 
x c (f) and controller inputs u c (t). Most linear time-invariant controllers can be put 
into this form with individual feedback gains appearing in the elements of the state 
matrices A c , B c , C c , and D c . Note that a closed-loop system is formed by connect- 
ing:^) the outputs y c (i) in equation (3.5) of the controller state model to the inputs 
u' p (t) in equations (3.1)-(3.3) of the plant model and, (2) the plant output vector y' a (t) 
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to the input vector u c (t) of the controller state model. This action of loop closure 
between the plant and the controller models is automatic and carried out inside the 
design package SANDY. 

In chapter 4 we will describe the formulation of a TECS controller using equations 
(3.4) and (3.5). 

3.2 Description of a Design Algorithm Based on Nonlinear Constrained Optimiza- 
tion 

As mentioned in the introduction, the design method proposed in this study is based 
entirely upon the minimization of a performance index for a set of closed-loop systems 
subjected possibly to additional design constraints such as closed-loop stability, co- 
variance responses to process and sensor noises, and bounds on robustness measures. 
A nonlinear programming technique [Ref.8] implemented in the numerical library 
NPSOL is used to solve the control problem with nonlinear design constraints. 

Typical performance measures in optimal control are quadratic penalties on the 
closed-loop system output and control responses to random disturbances. Three types 
of random disturbances are considered in the design problem: initial conditions, im- 
pulse inputs and white-noise inputs. Depending on the types of disturbances, different 
formulations of the performance index will be used. For the i th plant condition, the 
performance index to initial conditions or impulse inputs has the form 

M‘i) = 5 J‘‘ e 2 *' 1 E\y' p r (t)Q-yi(t) + u? ( l)R%{t)\dt (3.6) 

and for white-noise inputs the performance index is 

Wl) = \E..\x?(t,)Qy),(t l ) + 4 (3.7) 

where y l p (t) and txj,(f) are respectively the closed-loop responses of the performance 
and controller output variables (since u' p (t) = y c (t)) ■ The operator E ot [-] is the 
expected- value operator corresponding to a closed-loop system destabilized by a factor 
Qf t '. The matrices Q' and R' are the usual penalty weighting matrices on the system 
performance outputs and control inputs respectively. Evaluation of the performance 
index in equation (3.6) is done using the equivalence relation established in Appendix 



17 


A. The parameter a, is used to ensure that, when a steady-state optimal solution has 
been found, the closed-loop system eigenvalues of the controllable modes will have 
real parts less than —a, at the i th plant condition. The formulation of the outputs 
y' p {t) in equation (3.3) is general and can be defined to include the plant states xj,(t), 
the controller states x c (t) and the controller inputs u c (t). The plant states x' p (t) are 
modeled in the performance outputs y p (t) by letting C p = I, D p = fl p = 0. The 
controller states x c (t) and the controller inputs i i c (t) are modeled through the term 
D p in equation (3.3) and with u' p = y c established in the feedback connection. 

Design gains from selected parameters within the controller matrices are deter- 
mined such that a performance index of the form 

N p 

i=l 

is minimized. Individual performance index J,(z = 1 ,N P ) is one of those shown 
in equations (3.6) and (3.7). Performance indices are evaluated to an a-priori se- 
lected finite terminal time tj during the initial optimization. Efficient algorithms to 
evaluate the performance index J{tj ) and its gradients with respect to the design pa- 
rameters have been developed [Ref.6] that are comparable to those in the evaluation 
of the steady-state performance index (i.e based on Lyapunov solutions) in terms of 
computational time and accuracy. In fact, steady-state covariance responses can be 
accurately calculated by using a sufficiently large tj (e.g. at least 4-times the slowest 
time-constant of the closed-loop system modes). The terminal time tj is increased 
following each successful convergence in order to recover the steady-state solution. 
This gradual process of recovering the steady-state solution ensures that the final 
optimal solution will be a stabilizing one. The criterion for determining whether a 
steady-state solution has been reached is based on the condition that value of the 
optimized objective function will not vary by more than .1 percent from its previous 
value obtained for a smaller terminal time tj. Experiences gathered so far seem to 
justify the use of a performance index based on a finite terminal time. This procedure 
avoids the difficulties encountered by conventional methods based on the Lyapunov 
equation for the steady-state covariance responses, requiring special provision to han- 
dle the case when the controller design becomes destabilizing during the search. Note 
that by starting with a sufficiently small terminal time t j our design method can even 


(3.8) 
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be initiated with a controller design guess that is not stabilizing. 

The performance index J given in equation (3.8) enables designers to incorporate 
classical measures of design robustness or sensitivity to plant modeling uncertainties 
of the following types, 

1. Gain and phase margins in the control and sensor paths, 

2. Roll-off in the broken-loop frequency responses evaluated at the control and 
sensor feedback paths, 

3. Sensitivity of the closed-loop eigenvalues and the closed-loop performance mea- 
sures to perturbation in the plant model parameters. 

The first and second items are useful measures of robustness to unmodelled dynamics. 
The third item is a measure of robustness to uncertainties in the coefficients of the 
plant model. 

As seen in equation (3.8), our design algorithm for robust low-order controllers 
[Ref.6] involves an objective function that combines together several quadratic penal- 
ties of closed-loop system responses at different plant conditions. A simultaneous 
control-law synthesis to a set of multiple plant conditions allows the direct inclusion 
of robustness measures of type 3 into the design objective. Design to robustness 
measures of types 1 and 2 can be handled through the use of frequency-weighted 
performance functions and //°°-norm bounds on selected system transfer function 
matrices [Refs. 7, 19]. 

Specifically, robustness to plant parameter uncertainties or variation in design 
condition is achieved by using a set of quadratic performance indices that encompass, 
in addition to the one at the nominal design condition (i=l), other off-nominal design 
conditions weighted by the design parameters (z = 2,...,N P ) as in equation (3.8). 
To impose design robustness to unmodeled dynamics, one may define additional plant 
conditions that have performance indices representing penalties on control and sensor 
loop responses to high-pass noises injected separately into the control and sensor 
paths, or on the /f°°-norm of desired system transfer function matrices. 

Design for command following and tracking problems is done by formulating syn- 
thesis models and performance indices that penalize transient responses of the errors 



19 


between the actual and the commanded outputs. The commanded outputs (e.g fil- 
tered step, ramp, sinusoidal functions) are derived from a command generator model, 
a part of the augmented plant model, by using a combination of initial conditions 
and impulse inputs. If proportional and integral controls are desired in the feedback 
design, it is crucial that formulation of the design performance index reflects the need 
for an integral control action. The integral gains in the controller matrices tend to 
become ineffective (i.e converge to small values) when the objective function consists 
of responses to gaussian random disturbances of zero mean. In general, integral con- 
trol requirements must be defined in terms of a performance index that penalizes 
transient error responses of the controlled variables in the presence of parameterized 
random disturbance inputs where the generating time functions have nonzero con- 
stant steady-state components. That is, integral action is meaningful in the optimal 
design when it is used to compensate plant responses in the presence of parameterized 
random constant disturbances. 

Many types of parameterized random commands, such as impulse, step, and ramp 
inputs, are derived from responses of linear time-invariant shaping filters to random 
initial conditions or to impulse inputs of random magnitudes (Appendix B, section 
B.l). Outputs of these filters are interpreted as parameterized random commands 
whose magnitude y cm d is a random vector with zero mean E[y cm d ] = 0 and covariance 
E[ycmdyJmd\ = Ycmd- Elements of the magnitude vector y^d are uncorrelated when 
its covariance matrix Y cmd is diagonal. 

Random disturbances, such as Dryden or approximate Von Karman turbulences, 
are outputs of stable linear shaping filters excited by white-noise inputs w(t) (Ap- 
pendix B, section B.2) with zero mean .E[u;(<)] = 0 and covariance £’[ie(^)ie T (r)] = 
W p 6(t-r). 

Design requirements which cannot be easily handled through a quadratic perfor- 
mance index are enforced through linear constraints, nonlinear constraints, and direct 
bounds on the design parameters. Nonlinear constraints on the closed-loop stability, 
performance, and robustness are defined as follows, 

• Closed-loop stability of selected eigenvalues: 

Ci ^ Ci min » 1 , ..., fl) 

<7, < cr.nuuoti = l,...,n) 


( 3 . 9 ) 

( 3 . 10 ) 
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where £, is the damping ratio and <7, is the real part of the i th eigenvalue. For- 
mulation of these non-linear constraints is described in Appendix C. 

• Closed-loop mean-square responses: 

%*(</)) (3.11) 

where the parameters cr? . and define respectively the allowable lower 

pt min Pi max 1 J 

and upper bounds on the mean-square response of the i th plant output y Pt to 
random disturbances. 

• Desired //°°-norm bounds on closed-loop performance variables: 

K min 5: HGyi^O^) ||oo ^ K max (3.12) 

where ||G f yti; (jcj)|| 00 = sup w a {G yw (juj)) and (7 yt// (s) is the transfer function ma- 
trix of the closed-loop system between the plant disturbance inputs u?(s) and 
the plant outputs y p (s). 



Chapter 4 

TECS PROBLEM FORMULATION 


4.1 Formulation of Synthesis Models 


The NASA B-737 TSRV aircraft [Refs. 12,13,24] is used in this study of a total en- 
ergy control system design based on nonlinear constrained parameter optimization. 
Longitudinal aircraft models are generated for the landing-approach and cruise flight 
conditions described in tables E.l and E.2 (Appendix E). The linearized aircraft 
dynamic models have the form 


i p (t ) = A p x p (t) + B p u p (t) + T p w p (t) (4.1) 

y p (t) = C p x p (t) -|- D p u p (t) -f Q p w p (t) (4.2) 

where x p (t ) contains the aircraft rigid-body states, u p (t ) the elevator and throttle 
control inputs, w p (t) the wind components and the command inputs, and y p (t) the 
measurement and the performance criterion outputs. The aircraft models are subse- 
quently augmented with appropriate linear control actuation models of the form 


-^•oct(0 — A act x act (t) -J- B aC {U oc j(t) (4.3) 

U p(0 = (-'act^acti 0 T F) act ^acl(0 (4.4) 

Turbulence and command generation models are described by 


i<*(0 = AdXd(t) + B d w d (t ) (4.5) 

w p (t) = C d x d (t ) + D d w d (t) (4.6) 


They contain the commonly used Dryden or approximate Von Karman turbulence 
spectra [Ref.23]. Power spectral densities for Dryden turbulence models are given by 


$uM 


2 1 
Tu ^V l+(L u u/Vf 


2 F w 


1 + 3 (. L w u/vy 


'2 irV 


i + (L w u>/vy 


(4.7) 


^(w) = 


(4.8) 
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where <t u and a w are the root-mean-square intensities of turbulence (figure 4.2) with 


2 2 


L u 


(4.9) 


The parameters L u and L w are turbulence scale lengths defined as a function of 
altitude: 


• For altitude above 1750 ft: L u = L w = 1750 ft 


• For altitude below 1750 ft: L u = 145/t 1 ^ 3 , L w = h where h is the altitude in feet. 

Aircraft rigid-body dynamics are augmented with models for the control actuators 
and disturbances; the latter variables contain not only the wind inputs but also the 
command signals (Appendix B). The resulting plant model forms a synthesis model 
in the form of equations (3. 1 )— (3.3) . An automated procedure for formulating the 
synthesis model is implemented in a user-defined function for the MATLAB 1 control 
system software. A listing of the command procedure is given in Appendix F. 


4.2 An Optimal Design Approach for TECS 


The design algorithm for robust low-order controllers [Ref.6] implemented in the com- 
puter design package SANDY is applied to the design of an integrated AFCS. The 
total energy concept defines the “inner-loop” structure for an integrated autothrot- 
tle/autopilot control system and provides a basis for a multiloop control design. The 
basic TECS feedback controller shown in figure 2.2 with a pitch damper can be for- 
mulated into a controller state-space model of equations (3.4) and (3.5) as follows, 


1 iiE ^ 


" 0 0 0 ' 


1 XlE } 


' -1 -l/g 0 0 1 1/g ‘ 

X IL 

= 

0 0 0 


X IL 

+ 

1 -l/g 0 0-1 l/g 

iv c J 


0 0 0 


x v c ) 


0 0 0 0 0 1 


(4.10) 


1 MATLAB is a trademark of The Math Works, Inc. 
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UJ 


(4.11) 


Proportional and integral gains Kjp, Kpp, Kji and Kpi along with the pitch damper 
gains K q and I\$ and the additional gains k\ and &2 are design parameters in the 
controller state matrices. The gain schedule parameters Kgw and Kcas are assumed 
fixed at each plant condition and their contribution to the overall feedback gains can 
be absorbed into the other feedback gains. The design parameters are selected from 
the minimization of a performance index J of the form 

Np 

J = Y^W X J X (4.12) 

1=1 

where 

AT* 

•* = L-A“ , ('=l N r) < 413 * 

fc= 1 

and is a weighting factor assigned to the performance index J, at the i tk plant 
condition. Within each individual plant condition, performance indices (k — 
l,..., Nr) are used to address Nr different control design requirements. Design re- 
quirements defined in jj k ^ may be conflicting. Compromises among the conflicting 
requirements are usually made in the final design solution through iterative adjust- 
ment in the penalty weighting matrices. 

In the optimal TECS design problem, the objective function at each flight condition 

/ L \ 

is made up of at most three different performance indices J{ (k = 1,2,3), 


j<=iz j! 


k = 1 


(4.14) 
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The first performance index J} 1 is formulated to address design requirements in 
command tracking performance. It is given by 

•/!” = lim i f E[Q,(l(t) - 7 c(()) 2 + Oj(V( 0 - v;(t)) 2 \dt (4.15) 

tf-*oo l Jo 

This performance index is evaluated to parameterized random filtered step commands 
in V c and 7 C with V c {t) = V co ae~ at fi(t) and 7 C (0 = 7co(l — e~ at )fi(t) where ji{t) is the 
unit-step function. The parameter a determines the bandwidth of both the acceler- 
ation and the flight-path commands. In this design case, we use a = 1.2 rad/ sec, a 
typical value for flight- path and velocity command bandwidths. The variables V co and 
7 co are random parameters with zero means and covariances £[V(. 2 ] = &y = 1.0(/ps) 2 
and E[ 7 2 0 ] = <r 2 c = 1.0 (deg) 2 . The command V c (t) is computed simply as the integral 
of the acceleration command V c (t), i.e. 

v e (t)= fv c {r)dr (4.16) 

Jo 

The quantity is an output of the controller model given in equation (4.11). It 

is used in the criterion output V(t) — V*(t) and is defined as a linear combination of 
the commanded velocity V c (t) and flight-path angle 7 C (0> namely 

v;(t) = wt) + k 2lc (t) ( 4 . 17 ) 

The parameter k x depicts the achieved level of commanded velocity in steady-state. 
Recall that the TECS design structure in figure 2.2 does not include feedback of veloc- 
ity error, hence it is a type 0 system in the velocity variable V. Thus, in steady-state 
the aircraft velocity V(t) does not settle to the command value V c (t), and the param- 
eter k\ is always strictly less than one (i.e k\ < 1). The parameter k 2 is the amount of 
cross-coupling between flight-path command 7 c (t) and aircraft steady-state velocity 
V(t). Either parameter ki and/or k 2 can be set to lie within some desired values us- 
ing direct bound constraints. For our design, the inequality constraint — 1 < k 2 < 1 
ensures that for one degree of flight-path angle command 7 c (f), the change in air- 
craft steady-state velocity is less than one foot- per- second in magnitude. Clearly, 
for perfect decoupling where one does not allow changes in aircraft velocity during a 
flight-path angle command, k 2 would be constrained to be equal to zero. The parame- 
ters k\ and k 2 are determined during the optimization such that bound constraints on 
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k i and k^ are satisfied and the term (V(t) — V‘(t)) in the integrand of the performance 
index jj 1 * in equation (4.15) vanishes in the limit as tj — ► oo. 

The second performance index jj 2 ^ is set up to perform trade-off in the control 
bandwidth of the throttle loop. It is defined as 

Jf = lim i fWJJt,)] (4.18) 

tf->oo Z 

The performance index is evaluated to a high-pass noise input in the thrust 
command loop [Ref.25]. The noise input is generated from the response of a first- 
order filter to white-noise with zero mean £[£<^(2)] = 0 and covariance E[6 t ^ w {t + 
= The quantity 6 t k c (t) is the thrust feedback control as shown in 

figure 4.1. Cut-off frequency of the high-pass filter is set approximately equal to the 
desired broken-loop throttle control bandwidth (i.e. u>& thw = 0.2 rad/sec). Frequency- 
shaping of the disturbance input to the thrust command ensures that only control 
responses at high frequencies are penalized in j[ 2 \ 

Similarly, the performance index is used to perform trade-off in the control 
bandwidth of the elevator loop. It is given by 

y, (3 »=lin (4.19) 

In this case, J\ ’ is evaluated to a high-pass noise input in the elevator command loop. 
Again the noise input is obtained from the response of a first-order high-pass filter 
to white-noise with zero mean E[S ew (t)\ = 0 and covariance E[6 ew (t + T)<5 eu ,(t)] = 
a ew&( T )- The quantity £ ec (<) is the elevator feedback control. Cut-off frequency of 
the high-pass filter is approximately equal to the desired broken-loop elevator control 
bandwidth (i.e. us ew = 2.0 rad/sec). Note that the elevator control has a higher 
bandwidth than the throttle control. 

Other design considerations besides those depicted in the quadratic performance 
indices j[ l \ j[ 2 \ and jf 3 ^ are defined with the use of direct constraints. For example, 
desired closed-loop damping is achieved using the eigenvalue constraint defined in 
equation (C.41). Nonlinear damping constraints provide a direct means to achieve 
satisfactory damping of aircraft rigid-body modes. And covariance responses of se- 
lected outputs to clear air turbulence of Dryden spectra with RMS intensities from 
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the 99% probability level (figure 4.2) can be bounded using the covariance inequality 
constraint in equation (3.11). Nonlinear covariance constraints on control outputs 
ensure that the resulting optimal design has reasonable control activities to turbu- 
lence. 

4.3 Control System Design Tradeoff 

The design task involves a proper tradeoff among the following design performance 
goals: 

• Damping of system eigenvalues, 


• Command frequency response bandwidths, 


• Broken-loop crossover frequencies in the control paths, 


• Covariance responses of performance criteria, 


• Command tracking/following performance, 


• Disturbance rejection, 


• Robustness to plant parameter uncertainties and unmodeled dynamics. 

Performance objectives must be identified and established by the designer during the 
design tradeoff. If the desired objectives are not attainable, the NPSOL optimization 
[Ref. 8] would most likely fail to arrive at a feasible solution, e.g. design constraints 
are unattainable with the specified controller structure or design constraints are too 
stringent and need to be relaxed. 
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Minimum damping of the closed-loop system eigenvalues is obtained directly through 
specification of the nonlinear damping constraint £ > £ m ,„. If damping requirements 
are too stringent then performance in other areas will be sacrificed. We found that 
it is generally easier to achieve the desired damping level through the use of non- 
linear constraints on system modes than with output penalties in the cost function. 
Penalties on outputs which are partially related to lightly damped modes may un- 
necessarily degrade other system responses. Furthermore, outputs associated with 
lightly damped closed-loop modes will be different for different controller designs. 

Selection of desired command-loop frequency responses can be achieved through 
the penalties Q\ and Q 2 in the performance index jj 1 ^. Bandwidths of individual 
command-loop frequency responses K(s)/\£.(s) and l{s)/i c (s) relate to the ability 
of the system to track the respective commands. If the bandwidth in a particular 
command loop is too low, increasing the penalty Q{ ( i = 1,2) on the respective 
commanded output error in equation (4.15) will produce a higher bandwidth in that 
command path [Ref.25]. This systematic procedure enables designers to achieve sat- 
isfactory trade-off of performance in different command loops. 

Crossover frequency of a loop transfer function for a system broken at a control 
path defines approximately the control bandwidth of that path. Selection of the 
control-loop crossover frequencies is accomplished through the penalties R, ( i = 1,2) 
in equations (4.18) and (4.19). If the loop crossover frequency of a particular control 
path is too high, then higher penalty R , in the respective performance index, i.e. jj 2 ' 
or J{ , would lead to a design with a lower bandwidth [Ref.25]. Note that require- 
ments for simultaneous low control-loop crossover frequency and high command-loop 
bandwidth are generally contradictory. Usually, decreasing control-loop crossover 
frequencies will inevitably result in lower command-loop bandwidths. 

Upper bounds on root-mean-square (RMS) responses of selected performance cri- 
teria and control outputs to Dryden turbulence spectra may be specified to ensure 
adequate RMS responses. For the TECS design, throttle and elevator control activ- 
ities are the two variables of primary importance. Inequality bounds are therefore 
placed on covariance responses of 6 t h c {t) and 6 ec (t) to Dryden turbulence. The upper 
bounds are selected from the closed-loop control covariance responses of the current 
TECS design. In flight condition FLT2, one has to increase the elevator control activ- 
ity in turbulence by a small amount in order to achieve the desired level of closed-loop 
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damping. 

Design robustness in terms of closed-loop stability and performance measures to 
plant parameter uncertainties can be improved by using a performance index of the 
form given in equation (4.12), i.e. 

n p 

J = E W t J x 

i=i 

where each performance function J,(i = 2,...,Np) reflects design considerations at 
an off-nominal plant condition. Note that the index i = 1 corresponds to the nom- 
inal design condition. Stability in the presence of plant model parameter variations 
is defined in terms of damping constraints placed on closed-loop eigenvalues at the 
off-nominal plant conditions. When these constraints are satisfied, the resulting con- 
troller will then be robust to the given changes in plant conditions. However if the 
posed problem is overly constrained, then most likely a feasible solution cannot be 
found; in this case, a smaller set of uncertain plant models ought to be tried instead. 

An alternative formulation for the robustness criteria is through the //°°-norm 
bounds of selected system transfer function matrices, e.g ||G ? yw (ju;)|| 00 . The //°°-norm 
|| G '( ju ;)|| 00 is defined as the supremum of the maximum singular value of G(j u>) for 
u> 6 [0,oo). From the small gain theorem, H°° bounds on specific transfer function 
matrices give guaranteed gain and phase margins for nominally stable systems [Ref.3]. 
Nonlinear constraints of the form 


(4.20) 


||G yu ,(ju;)||oo <« ('1-21) 

can also be specified in the extended version of the SANDY design algorithm [Ref. 7] 
for robust stability and performance. Evaluation of this robust design procedure is 
left for a future study. 

4.4 Gain Schedule at Other Flight Conditions 

The design procedure described in section 4.3 will produce controller design gains 
optimal at one flight condition. Gain scheduling of the design at other flight conditions 
is usually required for optimum performance and may be achieved in one of three 


ways: 
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1. Design optimum controller gains at individual flight conditions separately. A 
gain schedule is developed for each controller gain. Usually it involves some form 
of curve-fitting with respect to airplane parameters such as calibrated airspeed 
or gross weight in the final design implementation. 

2. Preselect a gain schedule for each of the control or sensor paths. Design remain- 
ing controller parameters to satisfy all requirements with the preselected gain 
schedule structure over the entire set of design conditions simultaneously. 

3. Design a set of controller gains for a nominal flight condition. Redesign at 
other flight conditions based on gain scheduling a subset of controller gains. 
This procedure is simple and involves less computational effort in the actual 
implementation. 

The first method is straightforward and involves doing a separate design for each 
flight condition. However, this may result in a gain schedule which is too cumber- 
some to implement. In the TECS inner-loop structure, this method would require 
scheduling of six separate gains. 

The second approach involves deciding upon an a-priori gain schedule for the 
control or sensor paths, independent of the optimal design. These gains may be 
factored into the synthesis models at each flight condition. A single controller is 
then designed to satisfy all flight conditions with the pre-selected gain schedule. The 
design method involves defining a performance index similar to equation (4.20) that 
covers the entire range of flight conditions. It should be noted that this method 
depends upon the selection of a predefined gain schedule and provides no insight into 
how the gain schedule should be formulated. If the number of degrees-of-freedom 
is not adequately defined in the a priori selected controller structure, then most 
likely an optimal solution cannot be found that will meet all the design constraints. 
Thus, this method is generally not desireable since it sets the gain schedule around a 
predetermined structure that may be overly constrained. 

A more practical approach is to select a few gains in the controller which are to be 
gain scheduled. These gains are often introduced at the control input paths to provide 
compensation for changes in control effectiveness. The procedure in the third method 
begins with a set of gains optimized at a nominal flight condition perceived to be the 
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most important. Then the gain schedule is determined by optimizing only a selected 
subset of gains to other flight conditions while leaving the remaining controller gains 
fixed at their nominal values. For example, in the TECS inner-loop structure, the 
scheduled parameters are Kcas and Kgw while the gains Kep, Ktp, Kei, Kti, Ke, 
and I( q are fixed at the values optimized for the landing approach condition. The 
design procedure for optimizing the scheduled parameters and conducting design 
tradeoffs at each off-nominal flight condition is the same as described in section 4.3. 
This is a preferred approach since it provides insights into the level of improvement 
a gain schedule can offer with the selected degrees of freedom. 

4.5 Issues for Numerical Convergence 

Solution of the minimization problem is based upon a state-of-the-art nonlinear pro- 
gramming algorithm implemented in the NPSOL library [Ref.8]. Numerical condi- 
tioning of the TECS optimization problem is essential for successful convergence. 
Some key considerations are: 

• Convergence of the performance integral. The designers must verify that the in- 
tegrands of the quadratic performance index will approach zero in steady-state 
to the given set of input functions (i.e. white-noise, impulse, step, etc.). 


• Formulation of the constrained optimization problem. 


• Construction of the plant synthesis models. 


Performance indices in equations (3.6) and (3.7) and their gradients with respect 
to the design parameters are needed for optimization. The design algorithm imple- 
mented in the computer design package SANDY evaluates the performance indices 
and gradients to a finite terminal time tj. In steady-state (i.e. tj —* oo), the perfor- 
mance index in equation (3.6) or (3.7) would become unbounded when the closed-loop 
system contains neutrally stable or unstable modes that are both disturbable from 
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the input excitation and detectable in the performance function. Selection of an ap- 
propriately small finite terminal time tj would help avoid problems associated with 
lightly damped or unstable closed-loop poles in the initial phase of the numerical 
search. 

Design optimization to step commands or constant disturbance inputs requires 
careful formulation of the performance index, making sure that the integrands will 
settle to zero in the limit as tj — ► oo. Otherwise, the performance integral will be 
unbounded. For example, in the TECS problem formulation, the variable V* given 
in equation (4.17) was used instead of V c since the error (V — V c ) does not approach 
zero in steady state for a type 0 system and the performance integral of equation 
(4.15) would otherwise be unbounded. 

In problems dealing with constrained optimization, the possibility of having an 
overconstrained problem always exists. Careful formulation of the design objectives 
will definitely minimize problems encountered with nonlinear constraints and thereby 
speed up the design convergence. A helpful rule is to start the initial design with 
a minimal number of nonlinear constraints. Additional constraints can later be in- 
troduced in a systematic manner by order of relative importance and need. In the 
TECS design, damping in the aircraft rigid-body modes is of primary importance. 
Nonlinear constraints are henceforth used to achieve adequate closed-loop damping. 
Constraints on control covariance are included only when the initial design is found 
to be unsatisfactory. 

Modern nonlinear optimization techniques as found in the NPSOL library require 
that design parameters be of approximately the same order of magnitude. If some 
parameters are significantly larger than the others, then gradients of the cost function 
with respect to these parameters may appear too small; consequently the true opti- 
mum may not be reached to a given level of optimality tolerance. Proper selection 
of physical units in both the controller and the plant synthesis model ensures that, 
at the start, design parameters will have the same relative magnitudes. Here in our 
TECS problem, we found that the quantities 7 (rad) and V / g have the same order of 
magnitude and the feedback gains on these variables have the same size. 

Finally, it is crucial that the controller matrices in equations (3.4) and (3.5) form a 
minimal realization [Ref.6j. A nonminimal system will contain redundant parameters 
that lead to multiple local minima in the optimization. As a rule of thumb, the 
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maximum possible number n max of independent controller parameters is, for D c ^ 0, 

n m ax = (m + p)r + mp (4.22) 


and for D c = 0, 


nmax = (m + p)r 


(4.23) 


where r is the number of controller states, p is the number of controller inputs and m 
is the number of controller outputs. For example, in the TECS inner-loop structure as 
given in equation (4.11), a valid set of design parameters would be C c (l, 1), C c (2, 2), 
Z? c ( 1,1), D c ( 2,1), D c ( 2,3) and D c { 2,4) while the parameters D c { 1,2) and D c ( 2,2) 
are related to the above set through linear constraints, i.e Z) c (l,2) = D c (l, 1) and 
D c (2,2) = — Z) c (2, 1). 



Figure 4.1: High- Pass Frequency Shaping of Control Loop Activities 





Chapter 5 

INNER-LOOP DESIGN — SYNTHESIS AND ANALYSIS 


In this chapter we examine the synthesis of a TECS inner-loop design at two flight 
conditions FLTl and FLT2 (Appendix E). Results presented in the following sections 
illustrate the application of the design procedure described in section 4.2. 


5.1 Design for a Nominal Flight Condition 

For AFCS designs, tight airspeed and path control is important in the presence of 
windshear, particularly at low altitudes and during landing approach conditions. In 
this respect, the landing approach flight condition FLTl was considered most im- 
portant and is selected to be the nominal design condition for the TECS inner-loop 
gains. These TECS gains are designed to yield optimum performance at this flight 
condition. 

Designs at other flight conditions are synthesized around this nominal controller 
using the gain schedule parameters Kqas and Kqw (figure 2.2). When optimizing 
the TECS inner-loop gains, the gain schedule parameters are not considered as in- 
dependent parameters and must be fixed at some selected values. Without loss of 
generality, gain schedule at the flight condition FLTl is selected to be the same as 
that of the classical design. 

Design parameters in the TECS structure are determined following the procedure 
described in section 4.2. Arrival at a final design using constrained parameter op- 
timization involves compromise among different design performance goals. For the 
landing approach flight condition FLTl, the following design objectives were used: 


• The performance indices J[(i = 1,3) described in equations (4.15), (4.18) and 
(4.19) are given below. 

(a) Command frequency response bandwidths and command tracking perfor- 
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mance: 


= ,jim \ £ £[2(7(<) - iM? + (V(t) - v;(t)f)it (5.1) 

This performance index is evaluated to parameterized random filtered step com- 
mands in V c and 7 C with V c (t ) = V co ae~ at fi(t) and 7 C (<) = 7 C0 (1 -e _0 ‘)/i(<) where 
fj.(t) is the unit step function. The parameter a determines the bandwidth of 
both the acceleration and the flight-path commands. In this design case, we 
use a = 1.2 rad./ sec, a typical value for flight-path and velocity command band- 
widths. The variables V co and 7 co are random parameters with zero means and 
covariances E[V C 2 0 ] = a 2 Vc = 1.0 (fps) 2 and E[j 2 0 ] = a 2 c = 1.0 {deg) 2 . 

(b) Broken-loop crossover frequency in the throttle control loop: 

J \ 2) = t lim, ^£[ 10 *L(*/)] (5-2) 

The performance index is evaluated to a high-pass noise input in the thrust 
command loop. The noise input is generated from the response of a first-order 
high-pass filter to white-noise with zero mean = 0 and covariance 

a tkwH T )- I n this design, we choose a t h w = 1. The quantity 
MO ' s the thrust feedback control as shown in figure 4.1 . Cut-off frequency of 
the high-pass filter is set approximately equal to the desired broken-loop throttle 
control bandwidth (i.e. u!s lhv , = 0.2 rad/sec). 

(c) Broken-loop crossover frequency in the elevator control loop: 

J,' 3 ' = ^lhn \E[Z6l(t,)} (5.3) 

The performance index is evaluated to a high-pass noise input in the elevator 
command loop. The noise input is generated from the response of a first-order 
high-pass filter to white-noise with zero mean is[£ eu ,(<)] = 0 and covariance 
E[8 ew (t + T)6 ew (t)\ = In this design, we choose cr ew = 0.16 . The quantity 

8 ec (t) is the elevator feedback control as shown in figure 4.1 . Cut-off frequency of 
the high-pass filter is set approximately equal to the desired broken-loop elevator 
control bandwidth (i.e. = 2.0 rad/sec). 
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• Stability of the closed-loop eigenvalues: 

(a) Real part of the eigenvalues must be less than zero, 


<7, < 0 (i = 1, n) (5-4) 

(b) Damping ratio of the eigenvalues must be greater than 0.7, i.e 

c, >0.7(i = l,n) (5.5) 

• Mean square responses of control activities to clear air Dryden turbulence of 
cr u = 6.6 fps and a w = 6.3 fps, 

E[6 2 h J < 7.56 x I0~ 4 (lbst/lbsw) 2 (5.6) 

E[S] C ) < 6.735 x 10 -4 (rad) 2 (5.7) 

The above constraints will produce an optimal design (listed under SANDY 
in table 5.4) with a covariance in the speed variable V exceeding that of the 
classical design. To remedy this design problem, we include in the next design 
an additional constraint on the covariance response of the speed variable as 

E[V 2 } < 0.5122 {fps) 2 (5.8) 

The optimal design with this additional covariance constraint is shown in tables 
5.1 and 5.4 as the optimal design SANDY*. 

Sample of an input data file for the computer program SANDY corresponding to 
the above design formulation at flight condition FLT1 is given in Appendix F. The 
data file contains the plant synthesis model, controller model, disturbance specifica- 
tions, penalty weighting matrices, linear constraints, nonlinear constraints and direct 
bounds on parameters. The resulting optimal feedback gains and associated design 
parameters k\,k 2 are summarized in table 5.1. The following paragraphs discuss de- 
sign results obtained at flight condition FLT1 corresponding to the final optimized 
design SANDY*. 

In the optimal design SANDY*, the proportional gain Ktp to the throttle, pitch 
damper gain Kg and the decoupling parameter k 2 differ significantly from the classical 
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design. Improved performance and damping are achieved with a nonzero value of Ktp 
and a lower gain Kg. Recall that the magnitude of the parameter k 2 represents the 
degree of cross-coupling between flight-path command to steady-state speed error. 
This value is reduced in the optimal design in comparison to the classical design. 
Gains in the pitch damper have been reduced while those in the proportional and 
integral paths to the elevator are kept at the same level as the classical design; hence 
maintaining the desired command bandwidth and tracking performance. 

Table 5.2 gives the damping ratio and natural frequency of the closed-loop system 
eigenvalues. A minimum damping of 0.7 is achieved through the use of nonlinear 
constraints; this is an improvement over the minimum damping of 0.64 in the classical 
TECS design. 

Figures 5.1 and 5.2 show the command and broken-loop control frequency re- 
sponses of the classical and the optimized TECS designs respectively. Command and 
control-loop bandwidths for these designs are summarized in table 5.3. One design 
consideration is to have equal bandwidths in both the flight-path and acceleration 
command loops. Equal bandwidths in these command paths imply that transient 
errors in the total energy rate and energy distribution rate will decay exponentially 
to zero at a same time constant (i.e. coordinated flight-path and speed command re- 
sponses). The command bandwidths for the optimal design are similar to those of 
the classical design. Notice that the command bandwidths are achieved with a lower 
bandwidth in the elevator control path. 

Control-loop bandwidths give measures of control sensitivity to unmodeled high- 
frequency dynamics in the control paths. Lower control-loop bandwidth means less 
sensitivity to these types of unmodeled dynamics. The elevator loop bandwidth in the 
optimal design is smaller than that of the classical design at flight condition FLTl. 

Figures 5.3 to 5.6 show time responses of the closed-loop systems to step commands 
in flight- path angle and velocity. Velocity command is derived from an acceleration 
command as in equation (4.16). Results to a velocity command are shown in fig- 
ures 5.5 and 5.6. Both the flight-path angle and velocity commands are implemented 
through first-order command shaping filters with bandwidths equal to those used 
in the control-law synthesis (i.e. a = 1.2 rad/sec for both commands). Improved 
damping in the transient responses of flight-path angle and velocity errors is seen in 
figures 5.4 and 5.6. Responses to flight-path angle and velocity commands are similar 
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to step responses of a first-order filter. And there is less cross-coupling between the 
flight-path command and velocity. 

Table 5.4 compares the closed-loop rms responses of the classical and the optimal 
designs to Dryden clear air turbulence at the 99% probability level of intensities. 
Results clearly illustrate the effective usage of direct bounds on the respective output 
covariances. By imposing inequality constraints on the control covariances, one is 
guaranteed at the outset that the optimal design, when converged, will have the 
same or lower control activities than the classical design. Furthermore, covariance 
responses of other aircraft variables are lower or equal to those of the classical design. 
In particular, the covariance response of the speed variable V is equal to the value of 
the classical design due to the inequality constraint defined in equation (5.8). 

Conventional single-loop robustness 1 analysis gives allowable variations of gain or 
phase at either the control or sensor paths, while other remaining loops are closed 
at the nominal gains. Gain and phase margins determined from Bode or Nyquist 
plots yields the largest allowable variations (one-loop at a time) in loop gain and 
phase respectively starting from a nominally stable closed-loop system. Table 5.5 
summarizes values of single-loop stability margins. Single-loop stability margins are 
found to be satisfactory for both the optimal and classical designs. 

Multivariable robustness tests (Appendix D) based on singular values of loop 
return-difference matrices provide additional measures of design robustness in terms 
of guaranteed multivariable gain and phasfe margins [Ref.20]. The results may be con- 
servative but they are applicable to simultaneous variation of gain and phase in each 
loop. Figures 5.7-5.12 show plots of minimum singular values of the return-difference 
and the inverse-return-difference transfer matrices respectively at both the control 
and sensor paths. In these plots, diagonal scaling [Ref.23] on the transfer function 
matrices has been used to reduce conservatism and thereby improve estimates of ac- 
tual multiloop stability margins. Table 5.6 gives guaranteed multivariable stability 
margins in both the control actuator and sensor paths according to equations (D.8) 
and (D.14) of Appendix D. Both the optimal and classical designs possess similar 
multivariable stability margins. Note that requirements for robustness have not been 


1 Robustness margins are defined as the amount of allowable gain and phase increases prior to the 
onset of instability. 
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defined in the control-law synthesis. Systematic methods for improving design ro- 
bustness are discussed in section 4.2. Robust control-law synthesis using H°°-bounds 
or //-measure will be the subject of a future study. 

5.2 Determination of Gain Schedule 

Gain scheduling of the TECS inner-loops is necessary in order to satisfy performance 
and stability requirements at other flight conditions. Gain scheduling between a 
landing approach condition and one other cruise condition FLT2 is considered in this 
study. The objective is to demonstrate the usage of nonlinear constrained optimiza- 
tion in the design of gain schedule. Design results are discussed in the following 
paragraphs. The gain schedule design is obtained by re-optimizing the gains Kgw 
and Kqas in the throttle and elevator paths respectively. The TECS proportional, 
integral and pitch damper gains KTp,KEP,KTi,KEi,K q , and Kg are held constant 
at the values optimized for the landing- approach condition FLTl. Another set of 
design objectives are established for the cruise condition similar in form to that given 
in equations (5.1)-(5.8). They are as follows, 

• The performance indices J\{i = 1,3) described in equations (4.15), (4.18) and 
(4.19) are given below. 

(a) Command frequency response bandwidths and command tracking perfor- 
mance: 

Jf" - lim i /‘' £[2(7(i) - 7 C (0) 2 + 2(V(0 - V;(t))’]it (5.9) 

tf—oo l Jo 

This performance index is evaluated to parameterized random filtered step com- 
mands in V c and 7 C with V^(0 = K 0 ae -0 ‘/i(f) and 7 c (f) = 7 co (l — e - °‘)/i(f) 
where //(<) is the unit-step function. The parameter a determines the band- 
width of both the acceleration and the flight-path commands. In this design 
case, we use a = 1.2 rad/ sec, a typical value for flight-path and velocity com- 
mand bandwidths. The variables V co and 7^, are random parameters with zero 
means and covariances E[V^ = oy c = 1.0(/ps) 2 and E[ 7 2 0 ] = <r 2 e = 1.0 (deg) 2 . 

(b) Broken-loop crossover frequency in the throttle control loop: 

]?= lim i£[20<L(</)] 


(5.10) 
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/<j\ 

The performance index J{ ; is evaluated to a high-pass noise input in the thrust 
command loop. The noise input is generated from the response of a first-order 
high-pass filter to white-noise with zero mean £'[^ ( / lUJ (t)] = 0 and covariance 
+ = (T thwS( T )- Again, we choose a t hxu = 1. The quantity <W(<) 

is the thrust feedback control as shown in figure 4.1 . Cut-off frequency of the 
high-pass filter is set approximately equal to the desired broken-loop throttle 
control bandwidth (i.e. 0Js thw = 0.2 rad/sec). 

(c) Broken-loop crossover frequency in the elevator control loop: 

«A (3) = lim \E{U] c {tj)} (5.11) 

The performance index j| 3 ^ is evaluated to a high-pass noise input in the elevator 
command loop. The noise input is generated from the response of a first-order 
high-pass filter to white-noise with zero mean £’[<5 et „(t)] = 0 and covariance 
E[S ew (t + 7")<5 e u>(0] = a ew$( T )- Again, we choose a ew = 0.16 . The quantity 6 ec (t) 
is the elevator feedback control as shown in figure 4.1 . Cut-off frequency of the 
high-pass filter is set approximately equal to the desired broken- loop elevator 
control bandwidth (i.e. = 2.0 rad/sec). 

• Stability of the closed-loop eigenvalues: 

(a) Real part of the eigenvalues must be less than zero, 

<r , < 0 (i — 1, n) (5.12) 

(b) Damping ratio of the eigenvalues must be greater than 0.6, i.e 

(t > 0-6 (* = l,n) (5.13) 

• Mean square responses of control activities to clear air Dryden turbulence of 

= 4.7 fps and cr w = 4.7 fps , 

£[£ t 2 J < 1.804 x 10 ~ 5 {lbst/lbswf (5.14) 

E[S] C \ < 3.2 x 10 -3 (rad) 2 


(5.15) 
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It turns out that no feasible solution - can be found to satisfy simultaneously the 
closed-loop damping constraint of at least 0.6 and an elevator control activity to 
turbulence less than or equal to that achieved under the classical design. Thus 
the covariance bound on the elevator control is set to a value slightly higher than 
the value obtained from the classical design. Note that the selected bound on 
the elevator control activity to turbulence still yields acceptable level of control 
activity to the given level of turbulence. Recall that in the gain schedule design, 
the optimization has only two design parameters Kgw and Kcas ■ As we have 
demonstrated in flight condition FLTl, results of the classical design at flight 
condition FLT2 could similarly be achieved if we re-optimize all the inner-loop 
gains. Re-optimization of all the inner-loop gains at the flight condition FLT2 is 
presented in chapter 7 along with an altitude and speed-hold autopilot design. 

Gain schedule developed for the parameters Kgw and Kcas at the two design flight 
conditions and associated design parameters k\ and k? are shown in table 5.7. As 
previously mentioned, gain schedule at the nominal flight condition FLTl is chosen 
to be the same for both the optimal and classical designs. 

At flight condition FLT2 the optimized gain schedule based on the inner-loop 
design optimized at flight condition FLTl is quite different than the one chosen 
in the classical design. Design procedure for the selection of gain schedule in the 
classical design is usually based on closeddoop stability. It may not consider other 
design objectives such as tracking performance and control activities in the presence 
of turbulence. 

Table 5.8 shows the damping ratio and the natural frequency of the closed-loop 
system eigenvalues. A minimum damping of 0.6 is achieved in the optimal design 
through the use of nonlinear constraints on the closed-loop system eigenvalues. The 
result is improved over the conventional TECS closed-loop damping of 0.45. It is 
found by increasing the minimum damping requirement that higher damping of 0.7 
for all modes is not achievable with simple gain scheduling of the parameters Kgw 
and Kcas ■ Furthermore, the covariance constraint on the elevator control activity to 
turbulence with upper bound set equal to the value obtained in the classical design 
cannot be satisfied together with a minimal damping of 0.60 using only the design 
parameters Kgw and Kcas • 
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Figures 5.13 and 5.14 show the command and the broken-loop control frequency 
responses at flight condition FLT2. Command and control-loop bandwidths are sum- 
marized in table 5.9. Bandwidths in the flight-path and acceleration command loops 
are made nearly equal in the optimal design through adjustment of the design weight- 
ing parameters Qi and Q 2 in the cost function JqK The throttle control-loop band- 
width is significantly reduced in the optimal design; however, this result comes at 
the expense of an increased bandwidth in the elevator control path. The increase 
in elevator control bandwidth did not result in significantly higher elevator control 
activities to turbulence as seen in table 5.10. This design result illustrates the bene- 
fit of control covariance constraints in maintaining a given level of control activities 
while control bandwidths are optimized (i.e. increased) for improved damping and 
frequency of the longitudinal modes. 

Figures 5.15 to 5.18 show the time responses to step commands of the flight-path 
angle and velocity variables. Velocity command is developed from an acceleration 
command according to equation (4.16). Results associated with the velocity command 
are shown in figures 5.17 and 5.18. Both the flight-path angle and velocity commands 
are implemented through first-order command shaping filters with bandwidths equal 
to those used in the design synthesis (i.e. a = 1.2 rad/sec for both commands). 
Improved damping in the transient responses of flight-path angle and velocity errors 
is seen in figures 5.16 and 5.18. Command overshoots in flight-path angle responses 
have been nearly eliminated in the SANDY" design. 

Table 5.10 lists the closed-loop rms responses of the optimal and classical designs 
to Dryden clear air turbulence at the 99% probability level of intensities. Control rms 
responses are constrained with nonlinear inequality bounds on the control covariances. 
Results indicate that the optimal design possesses significantly lower throttle activi- 
ties. Nearly the same rms turbulence response in the elevator control is maintained 
in spite of significantly higher elevator control bandwidth. 

Table 5.11 gives values of single-loop type gain and phase margins. Single-loop 
stability margins are found to be satisfactory for both the optimal and the classical 
gain schedule designs. 

Figures 5.19-5.24 show respectively plots of minimum singular values of the return- 
difference and the inverse return-difference transfer matrices for the optimal and 
classical designs at both the control and sensor paths. Diagonal scaling has been 
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used to improve estimates of actual multivariable stability margins. Table 5.12 show 
guaranteed stability margins in both the control actuator and sensor paths obtained 
according to equations (D.8) and (D.14) in Appendix D. Note that the design robust- 
ness was achieved without direct design considerations. As pointed out in section 
4.2, H°° bounds on appropriate loop transfer matrices can be set up to address design 
robustness when such an improvement is needed. 

In the next chapters, we discuss the design of an altitude and speed-hold con- 
trol system using the TECS controller structure as inner-loop. The objective is to 
demonstrate the design features and the usage of nonlinear constrained parameter 
optimization in autopilot designs. There are two basic approaches one can adopt in 
the autopilot design: 

• The outer-loop design can be developed based on an existing (i.e previously 
designed) inner-loop TECS controller. This approach is simple and involves 
a fewer number of design parameters. But it does not yield the maximally 
achievable performance that would be possible if the inner-loop gains are re- 
optimized together with the outer-loop control functions. As seen in chapter 6, 
this method yields excessive throttle and elevator control activities to altitude 
command. 

• The outer-loop design is integrated with the synthesis of the inner-loop control- 
law. This latter procedure provides a better overall control design as demon- 
strated in chapter 7. 
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Table 5.1: Inner-Loop Feedback Gains (FLT1) 


Parameter 

SANDY 

SANDY' 

Classical 

I\tp 

0.09694 

0.14189 

0.0 

Kep 

3.3760 

3.4312 

3.36 

Kti 

0.3499 

0.3835 

0.4 

Kei 

2.4037 

2.3360 

2.52 


3.0071 

2.9788 

4.0 

Kb 

3.1240 

3.1184 

6.0 

k x I 

0.9194 

0.9060 

0.8932 

k 2 

-0.1633 

-0.0950 

0.2465 


Table 5.2: Closed-Loop System Poles (FLT1) 


Design 

SANDY ' 

Classical 

Mode 

c 

uj n (rad /sec) 

c 

u n (rad/sec) 

Phugoid mode 

0.7 

0.81 

0.64 

0.61 

Short period 

0.7 

2.23 

0.79 

3.01 


Table 5.3: Command and Broken-Loop Control Bandwidths (FLTl) 


Command/Control Path 

Bandwidths (rad/sec) 

SANDY' 

Classical 

7c 

0.8 

0.7 

Vc 

0.7 

0.8 

6" 

3.1 

4.3 

&thc 

0.37 

0.37 
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Table 5.4: Closed-Loop RMS Responses to Turbulence (FLT1) 


(<7 U = 6.6 fps, <j w = 6.3 fps) 


Variable 

SANDY 

SANDY' 

Classical 

7 (deg) 

0.9243 

0.9025 

1.0531 

V (fps) 

0.7915 

0.7157 

0.7157 j 

™zcg (§) 

0.0805 

0.0803 

0.0800 

6 t h (lbst/lbsw) 

0.02570 

0.02605 

0.02725 

6 e (deg) 

1.3324 

1.3311 

1.3324 


Table 5.5: Single-Loop Type Stability Margins (FLTl) 


Design 

SANDY' 

Classical 

Margins 

Gain Margin 

m 

Phase Margin 
(deg) 

Gain Margin 

(dB) 

Phase Margin 
(deg) 

Actuator Paths 

(-42, -boo) 

-60 

(-44, +31) 

-57 

Sensor Paths 

(-15, +14) 

(-64, +62) 

(-17, +20) 

(-63, +74) 


Table 5.6: Guaranteed Multivariable Stability Margins (FLTl) 


Design 

SANDY' 

Classical 

Margins 

Gain Margin 
(dB) 

Phase Margin 
(deg) 

Gain Margin 
(dB) 

Phase Margin 
(deg) 

Actuator Paths 

(-14.9, +5.2) 
(-5.3, +16.4) 

±50.2 

(-14.8, +5.2) 
(-5.1, +14.6) 

±48 

Sensor Paths 

(-2.8, +4.2) 

±22 

(-3.3, +5.4) 

±27 






























Table 5.7: Inner-Loop Gain Schedule and Design Parameters 


Parameter 

FLT1 

FLT2 

SANDY’ 

Classical 

SANDY’ 

Classical 

Kcas 

0.9322 

0.9322 

0.1950 

0.07006 

Kgw 

80,000 

80,000 

67,328 

80,000 ; 

ki 

0.9060 

0.8932 

0.9801 

0.9717 

h 

-0.0950 

0.2465 

-0.3208 

0.1795 


Table 5.8: Closed-Loop System Poles (FLT2) 


Design 

SANDY- 

Classical 

Mode 

C 

ui n (rad/sec) 

c 

u> n (rad/sec) 

Phugoid mode 

0.6 

0.64 

0.45 

0.42 

Short period 

0.68 

3.95 

0.47 

3.43 


Table 5.9: Command and Broken-Loop Control Bandwidths (FLT2) 


Command/Control Path 

Bandwidths (rad/sec) 

SANDY- 

Classical 

7c 

0.65 

0.55 

Vc 

0.65 

0.63 

S' 

4.3 

2.7 

StH 

0.29 

0.37 
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Table 5.10: Closed-Loop RMS Responses to Turbulence (FLT2) 


(<r u = 4.7 1 

'ps, a w = 4.7 fps) 

Variable 

SANDY' 

Classical 

7 (deg) 

0.2095 

0.2094 

V (fps) 

0.1369 

0.2060 

n zcg (S) 

0.0862 

0.0805 

6 t h (lbst/lbsw) 

0.00355 

0.004142 

6' (deg) 

0.2222 

0.1667 


Table 5.11: Single-Loop Type Stability Margins (FLT2) 


Design 

SANDY' 

Classical 

Margins 

Gain Margin 
(dB) 

Phase Margin 
(deg) 

Gain Margin 
(dB) 

Phase Margin 
(deg) 

Actuator Paths 

(-56, + 00 ) 

-83 

(-51, +30) 

-70 

Sensor Paths 

(-oo,21) 

(-68, +105) 

(-oo,24) 

(-48, +124) 


Table 5.12: Guaranteed Multivariable Stability Margins (FLT2) 


Design 

SANDY' 

Classical 

Margins 

Gain Margin 
(dB) 

Phase Margin 
(deg) 

Gain Margin 
(dB) 

Phase Margin 
(deg) 

Actuator Paths 

(-15, +5.2) 
(-5.4, +17.8) 

±51.6 

(-11.24, +4.7) 
(-5.4, +17.4) 

±51.2 1 

Sensor Paths 

(-3.6, 6.3) 

±29.9 

(-3.1, 5.0) 

±25.3 
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Figure 5.3: 7 C Command Responses of Classical Design (FLT1) 
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Figure 5.4: 7 C Command Responses of Optimal Design (FLT1) 
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Figure 5.9: Robustness at the Control Paths of Classical Design (FLT1) 
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Figure 5.10: Robustness at the Control Paths of Optimal Design (FLT1) 
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Figure 5.11: Robustness at the Sensor Paths of Classical Design (FLT1) 
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Figure 5.12: Robustness at the Sensor Paths of Optimal Design (FLTl) 
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Figure 5.15: 7 C Command Responses of Classical Design (FLT2) 
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Figure 5.16: 7 C Command Responses of Optimal Design (FLT2) 
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Figure 5.19: Robustness at the Control Paths of Classical Design (FLT2) 
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Figure 5.20: Robustness at the Control Paths of Optimal Design (FLT2) 
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Figure 5.21: Robustness at the Control Paths of Classical Design (FLT2) 
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Figure 5.22: Robustness at the Control Paths of Optimal Design (FLT2) 
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Figure 5.23: Robustness at the Sensor Paths of Classical Design (FLT2) 



Frequency (Rad/Sec) 


Figure 5.24: Robustness at the Sensor Paths of Optimal Design (FLT2) 





Chapter 6 

AIRSPEED AND ALTITUDE-HOLD AUTOPILOT 
AROUND AN OPTIMIZED TECS INNER LOOP 


In the previous chapter, we have described the use of a total energy concept in 
the design of an “inner-loop” for an integrated autothrottle/autopilot control system. 
This controller possesses essential and physically intuitive feedback paths to achieve 
the following design specifications : 

• Stability of closed-loop eigenvalues, 

• Minimum damping ratio of eigenvalues, 

• Acceptable mean square responses of control activities to turbulences, 

• Command frequency response bandwiths and command tracking performance. 

The following study covers the development of autopilot designs for an airspeed- 
and altitude-hold system using the TECS control-law as the inner loop. It was 
described in chapter 4 that the basic structure of the TECS controller does not include 
feedback of airspeed and/or altitude errors. The primary function of the inner-loop 
TECS (including the pitch damper) is to provide appropriate tracking bandwidths 
in the flight-path and longitudinal acceleration variables along with adequate closed- 
loop stability. To ensure that the airplane will indeed provide steady-state tracking 
of airspeed and altitude, we need to design an outer-loop that involves feedback 
of airspeed and altitude errors as shown in figure 6.1. In this chapter we present 
the synthesis of such an autopilot design around the existing optimal TECS inner 
loop discussed in chapter 5. In the next chapter 7, our design for the airspeed and 
altitude-hold system will also include the design parameters of the inner loop for 
improved tracking performance. This latter design approach for an integrated control- 
law synthesis covering simultaneously the inner and outer loops is clearly the preferred 
approach as indicated by the results of chapter 7. 
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In the airspeed-hold design, the airspeed error AV = V(<) — V c (t) (where I/(<) 
is the actual airspeed and V c (t ) is the commanded airspeed) is fedback through a 
feedback gain K v to the acceleration command V c of the TECS “inner loop” as shown 
in figure 6.1. While in the altitude-hold design, the altitude error Ah = h(t) — h c (t) 
(where h(t ) is the actual altitude and h c (t ) is the commanded altitude) is fedback 
through a feedback gain Kh to the flight-path command 7 C . The outer-loop feedback 
gains K v and K h are selected simultaneously from the minimization of the following 
performance index, 

J\= l™ 1 /"' E{Q„(V(t) - H(i)) 3 + <?»(*(<) - K(t)?]dt (6.1) 

subjected to an additional damping constraint of £ > Cmi„- The performance index J\ 
is evaluated to parameterized random filtered step commands K(f) = V co ( 1 — e~ at )fi(t) 
and h e (t) = h co ( 1 — e~ at )fi(t) where /i(t) is the unit-step function. The parameter a 
describes the bandwidth of both the airspeed and altitude commands. In this case, 
we use a = 1.2 rad/ sec similar to the bandwidths used for the inner-loop TECS design 
described in chapter 5. The variable V co is a random variable with zero mean and 
covariance E[V^ = ay c = l.O(fps) 2 . And the variable h co is a random variable with 
zero mean and covariance E\hl 0 ] = <r 2 c = 1.0(/<) 2 . 

Note that if we have only one design parameter, say K v in an airspeed-hold system 
or Kh in an altitude-hold system, then either the gain K v or Kh can be conveniently 
designed by simple root locus techniques. In fact, the root-locus method is very useful 
even in an optimization-based design setting. It provides us valuable information on 
the type of trade-offs occurred in the outer-loop design when the inner-loop is fixed 
at an existing controller design. In particular, one can examine by root locus the 
degradation of closed-loop damping when we add the outer loop. Root locus plots 
are shown in figures 6.2 and 6.3. As we observe from these plots, increasing either I\ v 
or Kh will eventually result in a reduction of the minimum achievable damping in the 
closed-loop system eigenvalues. Particularly the gain K v has always a destabilizing 
effect. Therefore one can no longer expect to retain the damping achieved by the 
inner-loop design when the outer- loop is closed. Trade-off between minimum damping 
of the closed-loop eigenvalues and command tracking performance is arrived at by 
varying the desired minimum damping £ m , n and the respective penalty weightings 
Q\ and Q 2 in the performance index Jj of equation (6.1). For example, we have to 
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reduce £ m i n from 0.7 to 0.6 in order to achieve satisfactory tracking performance. An 
outer-loop design is developed for each of the two flight conditions considered in this 
study. 

6.1 Outer-Loop Design at Flight Condition FLTl 

At flight condition FLTl, the performance index J\ used in the design of the outer- 
loop feedback gains K v and K h is given by 

J 3 = .Jjm \ j F £[(V(I) - V c (t)) 2 + 0.1 (/.(«) - K(t)f\dt (6.2) 

subjected to a damping constraint of ( > 0.6. The penalties Q v and Q h are deter- 
mined after a few design iterations and they are chosen to give comparable tracking 
performance between airspeed and altitude commands. For example, if transient re- 
sponses in the airspeed error is large then increasing Q v in equation 6.1 would improve 
the tracking of airspeed command. 

The final design gains are shown in Table 6.1. The following analysis have been 
carried out on the final design: 

• Frequency and damping ratio of closed-loop eigenvalues, 

• Time responses to airspeed and altitude commands, 

• Command and broken-loop control bandwidths, 

• Closed- loop rms responses to turbulences, 

• Single-loop type stability margins, 

• Multivariable stability margins. 

Notice that the inner-loop gains are fixed at the values obtained in chapter 5 for 
the optimal TECS controller at this flight condition. Table 6.2 gives the frequency 
and damping ratio of the closed-loop eigenvalues. Notice a slight reduction in the 
closed-loop damping as to be expected according to the root-locus plots of figure 6.2. 
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However in the next chapter, where the inner-loop gains are redesigned along with 
the outer-loop gains in an integrated manner, we are able to increase the damping 
back to at least 0.7. 

Figure 6.4 shows the command and broken-loop control frequency responses while 
table 6.3 summarizes the command and broken-loop control bandwidths. The control- 
loop bandwidths are essentially the same as those obtained in the inner-loop design. 
However, as expected, the command bandwidths of the outer loop is generally lower 
than those of the inner loop. 

Table 6.4 gives the rms responses of the closed-loop system to Dryden clear air 
turbulence at the 99% probability level of intensities. Results are comparable to those 
achieved by the optimized TECS inner-loop design. 

Table 6.5 gives the single-loop type gain and phase margins at the actuator and 
sensor paths. The corresponding multivariable stability margins are shown in ta- 
ble 6.6. All these robustness results are found to be satisfactory. 

Figures 6.5 and 6.6 show the time simulation to step commands in altitude of 
1000// and airspeed of 10 fps respectively. Note that in both cases the commanded 
value has been reached in about 100 seconds. Evaluation of time responses to the 
respective commands provide insights into the iterative design procedure since these 
responses reflect directly the effects of varying the penalties Q v and Qh- It is seen 
that comparable performance is achieved between speed and altitude commands. One 
advantage of our outer-loop design is that there is small coupling between the speed 
and altitude variables. That is, there is little change in altitude when we apply an 
airspeed command and vice versa. Both the elevator and throttle control excursions 
in the altitude command responses are rather excessive. This is due to the fact that 
we have achieved a tight tracking on the altitude response and the outer loop has 
been synthesized around an existing (i.e fixed) inner loop. 


6.2 Outer- Loop Design at Flight Condition FLT2 

The discussion of the outer-loop design at the cruise condition FLT2 follows an outline 
similar to the one presented for the landing approach condition FLTl. At flight con- 
dition FLT2, the performance index J\ used in the design of the outer-loop feedback 
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gains K v and Kh is given by 

h = Jim l - j E[\0(V(t) - V c (t )) 2 + (h(t) - ft c (0) 2 ]d< (6.3) 

subjected to a damping constraint of ( > 0.6. 

The final design gains are shown in Table 6.7. The following analysis have also 
been carried out on the final design: 

• Frequency and damping ratio of closed-loop eigenvalues, 

• Time responses to airspeed and altitude commands, 

• Command and broken- loop control bandwidths, 

• Closed-loop rms responses to turbulences, 

• Single-loop type stability margins, 

• Multivariable stability margins. 

Notice that the inner-loop gains are fixed at the values obtained in chapter 5 for the 
optimal TECS controller at this flight condition. Table 6.8 gives the frequency and 
damping ratio of the closed-loop eigenvalues. In this design, we are able to maintain 
the same level of closed-loop damping of 0.6 as the inner-loop gain schedule design 
at the flight condition FLT2. However, in the chapter 7 where we consider the re- 
design of the inner-loop gains together with the outer-loop gains, an improvement in 
closed-loop damping of 0.65 can be obtained. 

Figure 6.10 shows the command and broken-loop control frequency responses while 
table 6.9 summarizes the command and broken-loop control bandwidths. The control- 
loop bandwidths are essentially the same as those obtained in the inner-loop design. 
However, as expected, the command bandwidths of the outer loop are generally lower 
than those of the inner loop. 

Table 6.10 gives the rms responses of the closed-loop system to Dryden clear air 
turbulence at the 99% probability level of intensities. Results are similar to those 
achieved by the TECS inner-loop design. 
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Table 6.11 gives the single-loop type gain and phase margins at the actuator and 
sensor paths. The corresponding multivariable stability margins are shown in ta- 
ble 6.12. All these robustness results are found to be satisfactory. 

Figures 6.11 and 6.12 show the time simulation to step commands in altitude of 
1000/i and airspeed of 10 fps respectively. Note that in both cases the commanded 
value has been reached in about 100 seconds. Evaluation of time responses to the 
respective commands provide insights into the iterative design procedure since these 
responses reflect directly the effects of varying the penalties Q v and Qh- It is seen 
that comparable performance is achieved between speed and altitude commands. One 
advantage of our outer-loop design is that there is small coupling between the speed 
and altitude variables. That is, there is little change in altitude when we apply an 
airspeed command and vice versa. 

In this chapter we have developed a design procedure applied to the synthesis of 
an autopilot design using an a-priori designed inner loop. A performance index based 
simply on tracking errors to step commands in airspeed and altitude was used. We 
are able to achieve good tracking performance and adequate closed- loop damping 
requirements. In some situation, it may be required to re-design the inner loop at 
the same time one is designing the outer loop in order to accomodate more stringent 
requirements in closed-loop damping and control activities to turbulence. The next 
chapter, we consider the problem of simultaneous inner and outer loop design in an 
airspeed and altitude hold autopilot. 
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Table 6.1: Outer-Loop Gains with Existing Optimal TECS Design (FLTl) 


Parameter 

SANDY 

Ktp 

0.1419 

Kep 

3.4312 

Kti 

0.3835 

Kei 

2.3360 


2.9788 

K b 

3.1184 

K v 

-7.5202 xlO -3 

I<h 

-1.1223 xl0~ 3 

Kcas 

0.9322 

Kgw 

80,000 


Table 6.2: Closed-Loop System Poles (FLTl) 


Design 

SANDY I 

Mode 

C 

u n (rad/sec) 

Phugoid mode 

0.6 

0.43 

Short period 

0.7 

2.23 

Altitude Mode 

1 

0.124 

Integral Speed Mode 

0.6 

0.37 


Table 6.3: Command and Broken- Loop Control Bandwidths (FLTl) 


Command/ Control Path 

SANDY 


Bandwidths (rad/sec) 

h c 

0.6 

V e 

0.35 

See 

3.19 

Sthc 

0.40 




Table 6.4: Closed-Loop RMS Responses to Turbulence (FLTl) 


(gu = 6.6 fps, 0 + = 6.3 fps) 


Variable 

SANDY 

7 (deg) 

1.2491 

V (fps) 

0.5487 

h (ft) 

11.3 

n zcg (g) 

0.09335 

6th (lbst/lbsw) 

0.03558 

6 e (deg) 

1.6901 


Table 6.5: Single- Loop Type Stability Margins (FLTl) 


Design 

SANDY 

Margins 

Gain Margin 
(dB) 

Phase Margin 
(deg) 

Actuator Paths 

(—19.6, Too) 

-42 

Sensor Paths 

(-8.3, +8.6) 

(-43, +62) 


Table 6.6: Guaranteed Multivariable Stability Margins (FLTl) 


Design 

SANDY 

Margins 

Gain Margin 
(dB) 

Phase Margin 
(deg) 

Actuator Paths 

(-8.5, +4.2) 
(-4.4, +9.3) 

±38.5 

Sensor Paths 

(-2.2, +2.9) 

±16.4 
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Table 6.7: Outer-Loop Gains with Existing Optimal TECS Design (FLT2) 


Parameter 

SANDY 

K T p 

0.10895 

Kep 

9.8268 

Kti 

0.29446 

Kei 

6.6903 

K, 

8.5310 

Ke 

8.9311 

Kv 

-2.1231 xlO -3 

I<h 

-2.1016 xlO -4 

Kcas 

0.07006 

Kgw 

80,000 


Table 6.8: Closed-Loop System Poles (FLT2) 


Design 

SANDY I 

Mode 

C 

(rad /sec) 

Phugoid mode 

0.60 

0.508 

Short period 

0.68 

3.95 

Altitude Mode 

1 

0.08578 

Integral Speed Mode 

0.8 

0.26 


Table 6.9: Command and Broken-Loop Control Bandwidths (FLT2) 


Command/Control Path 

SANDY 


Bandwidths (rad/sec) 

K 

0.10 

v c 

0.30 

6'c 

4.26 

&thc 

0.29 





Table 6.10: Closed-Loop RMS Responses to Turbulence (FLT2) 


(gu = 4.7 fps, <T W = 4.7 fps) 


Variable 

SANDY 

7 (deg) 

0.2596 

V (fps) 

0.1438 

h (ft) 

8.1443 

n zcg (g) 

0.09089 

6 t h (lbst/lbsw) 

0.004549 

<5 e (deg) 

0.2418 


Table 6.11: Single-Loop Type Stability Margins (FLT2) 


Design 

SANDY 

Margins 

Gain Margin 
(dB) 

Phase Margin 
(deg) 

Actuator Paths 

(-16, Too) 

-58 

Sensor Paths 

(-12, +11.4) 

(-50, +81) 


Table 6.12: Guaranteed Multivariable Stability Margins (FLT2) 


Design 

SANDY 

Margins 

Gain Margin 
(dB) 

Phase Margin 
(deg) 

Actuator Paths 

(-11, +4.6) 
(-5.1, +14) 

±47 

Sensor Paths 

(-2.6, +3.6) 

±20 
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Figure 6.1: Airspeed and Altitude-Hold Design with an Optimized TECS Inner Loop 
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Real Part (rad/sec) 

Figure 6.3: Root-Locus Plots versus Design Gains K v and A'/, (FLT2) 
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Figure 6.4: Loop Frequency Responses of Optimal Outer-Loop Design (FLTl) 















76 



Figure 6.7: Robustness at the Control Paths of Optimal Outer-Loop Design (FLTl) 



Figure 6.8: Robustness at the Control Paths of Optimal Outer-Loop Design (FLTl) 
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Figure 6.9: Robustness at the Sensor Paths of Optimal Outer- Loop Design (FLTl) 



Figure 6.10: Loop Frequency Responses of Optimal Outer-Loop Design (FLT2) 
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Figure 6.11: h c Command Responses of Optimal Outer-Loop Design (FLT2) 
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Chapter 7 

SIMULTANEOUS DESIGN OF TECS INNER AND 

OUTER LOOPS 


In chapter 5 we have described a design procedure for the synthesis of a TECS 
controller at two flight conditions FLTl and FLT2. In these inner-loop designs, con- 
troller design parameters were determined so that the closed-loop systems satisfy the 
desired stability, bandwidth and turbulence requirements. In chapter 6, an airspeed 
and altitude-hold autopilot design was considered using the existing optimized TECS 
inner-loop design parameters. In this chapter, we re-examine the autopilot design 
from the point of view of improving its performance beyond that achieved by the 
outer- loop design of chapter 6. 

Some of the areas that one can improve by re-designing the TECS inner loop are: 

• Closed- loop damping of system eigenvalues, 

• Closed-loop rms responses to turbulence, 

• Control activities to commands, 

• Design robustness. 

The control design structure consists of the TECS control-laws augmented with, in 
the outer loop, feedbacks of velocity and altitude errors as shown in figure 6.1. Note 
that we adopt the design philosophy that errors in airspeed will be fedback to the 
acceleration command V c through a feedback gain K v and altitude errors are corrected 
through feedback signals to the flight-path command 7 C via the gain K^. A state- 
space representation of this controller is given below: 
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(7.2) 


where AV^ = V — V c is the velocity error in fps and Ah = h — h c is the altitude 
error in ft. Design parameters in the simultaneous inner and outer loop design are 
the proportional and integral gains I\ T p, Kti , K E p and Ae/. I* 16 P ltc h damper gams 
K q and K b , the velocity-error gain I< v and the altitude-error gain K h . For each 
flight condition we will determine a new set of feedback gains, hence without loss of 
generality we use the same gain schedule values for K G w and I< C as in the control-law 
synthesis. In fact, their contribution to the overall feedback gains are accounted for 
by the other feedback gains. 

Following the design procedure described in chapter 4, the objective function for 
the airspeed and altitude-hold design consists of the sum of three performance indices 

#(*= 1 , 2 , 3 ), 

3 

J. = £ Jt 

k= 1 


(7.3) 
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The first performance index is formulated to address design requirements in 
airspeed and altitude command tracking performance. It is given by 

J i (1) = J“22o \ C E ^ v ^) ~ V 'M 2 + Q 2 W) - h c (t)) 2 }dt (7.4) 

This performance index is evaluated to parameterized random filtered step commands 
in V c and h c with V c (t) = V co (l - e~ at )fi{t) and h c (t) = h co (l - e where p(t) 
is the unit-step function. The parameter a determines the bandwidth of both the 
airspeed and the altitude commands. In this design case, we use a = 1.2 rad/ sec, 
the same value as the flight-path and velocity command bandwidths. The variables 
V co and h co are random parameters with zero means and covariances E[V 2 ] = = 

1 .0( Jpsy and £[/£] = al = 1 .0 (/()’. 

Note that this autopilot design formulation includes feedback of velocity error, 
hence making the system type 1 in the velocity variable V. In steady-state the aircraft 
velocity V(t) will settle to the command value V c (t) if the feedback gain I\ v is nonzero 
and stabilizing. Similarly, the altitude variable introduce an additional integrator 
in the altitude loop hence making it also a type 1 system. The altitude response 
will therefore settle to the commanded value if the system has a nonzero stabilizing 
feedback gain A*. Thus the integrand of the performance index in equation (7.4) 
vanishes in the limit as tj — > oo. 

The second performance index is set up to perform trade-off in the control 
bandwidth of the throttle loop. It is defined as 

J f 2) = (7.5) 

The performance index jJ 2 * is evaluated to a high-pass noise input in the thrust 
command loop. The noise input is generated from the response of a first-order high- 
pass filter to white-noise with zero mean £ l [<5 < ^ u ,(/)] = 0 and covariance E[S ( h w (t + 
r )^<Aui(0] ^thw^{ T )- The quantity ^tj, c (t) is the thrust feedback control as shown in 

figure 4.1. Cut-off frequency of the high-pass filter is set approximately equal to the 
desired broken-loop throttle control bandwidth (i.e. = 0.2 rad/sec). Frequency- 

shaping of the disturbance input to the thrust command ensures that only control 
responses at high frequencies are penalized in j[ 2 \ 
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Similarly, the performance index Jj 3 ^ is used to perform trade-off in the control 
bandwidth of the elevator loop. It is given by 


J| 3) = 


\\m l -E[R 2 8] c {tj)\ 

tf->oo l 


(7.6) 


. /o\ 

In this case, is evaluated to a high-pass noise input in the elevator command loop. 
Again the noise input is obtained from the response of a first-order high-pass filter 
to white-noise with zero mean E[8 ew {t)] = 0 and covariance E[8 ew (t + r)<5 eu/ (<)] = 
a ew^{ T )- The quantity 8 fC (t) is the elevator feedback control. Cut-off frequency of 
the high-pass filter is approximately equal to the desired broken-loop elevator control 
bandwidth (i.e. = 2.0 rad/sec at flight condition FLTl and u> Sew = 0.1 rad/sec at 

flight condition FLT2). The lower bandwidth selected at the cruise condition FLT2 
ensures that the elevator control bandwidth is of the same size as previous inner-loop 
design. 

Other design considerations besides those depicted in the quadratic performance 
indices ./{ l \ jj 2 \ and Jj 3 ' are defined as before with the use of direct constraints. 
For example, desired closed-loop damping is achieved using the eigenvalue constraint 
defined in equation (C.41). Nonlinear damping constraints provide a direct means to 
achieve satisfactory damping of aircraft rigid-body modes. And covariance responses 
of selected outputs to clear air turbulence of Dryden spectra with RMS intensities 
from the 99% probability level (figure 4.2) can be bounded using the covariance 
inequality constraint in equation (3.11). Nonlinear covariance constraints on control 
outputs ensure that the resulting optimal design has reasonable control activities to 
turbulence. 

In this design case, selection of initial guess for the design parameters and the 
penalty weighting matrices is facilitated by the results obtained in chapter 5 for 
the TECS inner-loop designs. All design gains were initially selected to be those 
determined in chapters 5 and 6 since these solutions have shown to provide reasonable 
stabilizing designs. Weighting matrices are also initially chosen to be those used in 
the previous design cases. Starting from these, one can iterate until a satisfactory 
autopilot design is attained. 
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7.1 Simultaneous Design of TECS Inner and Outer Loops at Flight Condition FLTl 

For the landing approach condition FLTl, the following design objectives were used: 

• The performance indices J\{i = 1,3) described in equations (7.4), (7.5) and (7.6) 
are given below. 

(a) Command frequency response bandwidths and command tracking perfor- 
mance: 

j! ‘ l= \ C £ [ ,0 ( V <‘) - + (*(<)“ MO) 5 ]* (7.7) 

This performance index is evaluated to parameterized random filtered step com- 
mands in V c and h c with V c (t) = V C0 {1 - e~ ai )p{t) and h c (t) = h co { 1 - e~ at )p(t) 
where p(t) is the unit-step function. The parameter a determines the band- 
width of both the airspeed and the altitude commands. In this design case, we 
use a 1.2 rad/ sec, the same value as the flight-path and velocity command 

bandwidths. The variables V co and h co are random parameters with zero means 
and covariances E[V C 2 0 ] = = 1.0(/ps) 2 and E[h 2 J = o 2 c = 1.0 (ft) 2 . 

(b) Broken-loop crossover frequency in the throttle control loop: 

J i (2) = Jjm, (7.8) 

The performance index is evaluated to a high-pass noise input in the thrust 
command loop. The noise input is generated from the response of a first-order 
high-pass filter to white-noise with zero mean £[<5t/iu,(f)] = 0 and covariance 
E[S lhw (t + r)6 (/lu ,(f)] = cr 2 hw 6(r). In this design, we choose (T thw = 1. The quantity 
&thc(t) is the thrust feedback control as shown in figure 4.1 . Cut-off frequency of 
the high-pass filter is set approximately equal to the desired broken-loop throttle 
control bandwidth (i.e. u s , hvi = 0.2 rad/sec). 

(c) Broken-loop crossover frequency in the elevator control loop: 

- 7 '' 31 = f £('/)] ( 7 . 9 ) 

The performance index is evaluated to a high-pass noise input in the elevator 
command loop. The noise input is generated from the response of a first-order 
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high-pass filter to white-noise with zero mean £[<5 etl ,(<)] = 0 and covariance 
•£'[^eu;(^ + r )<5 e w(0] = cr eu-^( 7 )- I n this design, we choose a ew = 0.16 . The quantity 
£ ec (0 is the elevator feedback control as shown in figure 4.1 . Cut-off frequency of 
the high-pass filter is set approximately equal to the desired broken-loop elevator 
control bandwidth (i.e. = 2.0 rad/sec). 

• Stability of the closed-loop eigenvalues: 

(a) Real part of the eigenvalues must be less than zero, 


a , < 0 (z = 1, n) 


(7.10) 


(b) Damping ratio of the eigenvalues must be greater than 0.7, i.e 


C, > 0.7 (i = l,n) 


(7.11) 


• Mean square responses of control activities to clear air Dryden turbulence of 
a u = 6.6 fps and a w = 6.3 fps, 

£’[<5^ c ] < 7.56 x \0~ 4 (Ibst/lbsw) 2 (712) 

E[8 ] J < 6.735 x 10" 4 (rad) 2 (7.13) 

The resulting set of optimal feedback gains is shown in table 7.1. Overall these gains 
differ significantly from those where the inner and outer loops are designed separately. 
Improved performance and damping are achieved with higher values of Kjp, I\ q , Ko 
and K v . Table 7.2 gives the damping ratio and natural frequency of the closed- 
loop system eigenvalues. A minimum damping of 0.7 is achieved through the use of 
nonlinear constraints; this is an improvement over the minimum damping of 0.60 in 
the optimal outer-loop design in chapter 6. 

Figures 7.1 shows the command and broken-loop control frequency responses of the 
optimized design. Command and control-loop bandwidths for these designs are sum- 
marized in table 7.3. Note that the design has slightly higher command bandwidth 
in the velocity variable than the outer-loop design in chapter 6, while the altitude 
command bandwidth is significantly reduced. This is primarily due to the fact that 
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we have limited the control activities to the values defined by the classical TECS 
design described in chapter 5 and we are imposing a higher damping requirement. 
Control-loop bandwidths are similar in the two designs. 

Figures 7.2 and 7.3 show time responses of the closed-loop systems to a step com- 
mand in altitude of 1000// and in airspeed of 10/ps respectively. Control activities 
are less than half of those shown in table 6.5 for the previous outer-loop design. Ta- 
ble 7.4 summarizes the closed-loop rms responses to Dryden clear air turbulence at 
the 99% probability level of intensities. Results clearly illustrate the effective usage 
of direct bounds on the respective output covariances. By imposing inequality con- 
straints on the control covariances, one is guaranteed at the outset that the optimal 
design, when converged, will have the same or lower control activities than the clas- 
sical design. Furthermore, covariance responses of other aircraft variables are lower 
or equal to those of the classical design shown in table 5.4. 

Table 7.5 summarizes values of single-loop stability margins. Single-loop stability 
margins are found to be satisfactory and improved over those in table 6.5 of the 
previously optimized outer-loop design. 

Figures 7. 4-7. 6 show plots of minimum singular values of the return-difference and 
the inverse-return-difference transfer matrices respectively at both the control and 
sensor paths. In these plots, diagonal scaling on the transfer function matrices has 
been used to reduce conservatism and thereby improve estimates of actual multiloop 
stability margins. Table 7.6 gives guaranteed multivariable stability margins in both 
the control actuator and sensor paths according to equations (D.8) and (D.14) of 
Appendix D. Again these results are better than those in table 6.6 for the optimized 
outer-loop design of chapter 6. 

7.2 Simultaneous Design of TECS Inner and Outer Loops at Flight Condition FLT2 

Instead of performing a simple gain schedule design based on the results obtained in 
the previous section 7.1 for the landing approach condition, we choose to re-design 
all the TECS inner and outer-loop gains at the cruise condition FLT2. The objective 
is to see whether one can improve the design performance beyond those achieved in 
chapter 5 and 6 combined. Another set of design objectives are established for the 
cruise condition similar in form to that given in equations (7.7)-(7.13). They are as 
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follows, 

• The performance indices J[{i = 1,3) described in equations (7.4), (7.5) and (7.6) 
are given below. 

(a) Command frequency response bandwidths and command tracking perfor- 
mance: 


J ’" = mv(l) - v <(‘)f + 2(h(t) - K(t)f)dt (7.14) 

This performance index is evaluated to parameterized random filtered step com- 
mands in V c and h e with V e (t) = V co (l - e~ at )fi{t) and h c {t) = h co (l - e~ at ) f i(t) 
where //(<) is the unit-step function. The parameter a determines the band- 
width of both the airspeed and the altitude commands. In this design case, we 
use a = 1.2 rad/ sec, the same value as the flight-path and velocity command 
bandwidths. The variables V co and h co are random parameters with zero means 
and covariances E[V c l } = a\ c = 1.0(/ps) 2 and E[h 2 J = a\ c = 1.0(/<) 2 . 

(b) Broken-loop crossover frequency in the throttle control loop: 

J,(2) = tjil^l E l 0 - 03S L((/)l (7.15) 

The performance index jj 2 * is evaluated to a high-pass noise input in the thrust 
command loop. The noise input is generated from the response of a first-order 
high-pass filter to white-noise with zero mean £[6 tW (0] = 0 and covariance 
£’[<W(^ + 'r)<^/ ltl ,(0] = ^Jhw^( T )- 1° this design, we choose a thw = 1. The quantity 
<W(0 ' s the thrust feedback control as shown in figure 4.1 . Cut-off frequency of 
the high-pass filter is set approximately equal to the desired broken-loop throttle 
control bandwidth (i.e. u}g ttiui = 0.2 rad/sec). 

(c) Broken-loop crossover frequency in the elevator control loop: 

J ' 3) = ^I 10(W ec(*/)] (7.16) 

The performance index jf 3 ^ is evaluated to a high-pass noise input in the elevator 
command loop. The noise input is generated from the response of a first-order 
high-pass filter to white-noise with zero mean £[6 eur (t)J = 0 and covariance 
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E[S ew (t + r)6 etl/ (f)] = t r* w 6(r ). In this design, we choose a ew = 1.0 . The quantity 
S ec (t) is the elevator feedback control as shown in figure 4.1 . Cut-off frequency of 
the high-pass filter is set approximately equal to the desired broken-loop elevator 
control bandwidth (i.e. u>$ tw =0.1 rad/sec). 

• Stability of the closed-loop eigenvalues: 

(a) Real part of the eigenvalues must be less than zero, 


<r, < 0 (* = l,n) (7.17) 

(b) Damping ratio of the eigenvalues must be greater than 0.65, i.e 

C, > 0.65 (i = l,n) (7. IS) 

Note that no feasible solutions can be found that will satisfy damping require- 
ment exceeding 0.65 and at the same time maintaining the covariance responses 
achieved under the classical design (Table 5.10). 

• Mean square responses of control activities to clear air Dryden turbulence of 
= 4.7 fps and a w = 4.7 fps, 

We are restricting our design to provide the same levels of covariance responses 
as the classical design shown in table 5.10. 


£[6« 2 J < 0.18037 x 10~ 4 {Ibst/lbsw) 2 

(7.19) 

E[S 2 J < 0.17581 x 10" 2 (rad) 2 

(7.20) 

E[V 2 } < 0.0424(/ps) 2 

(7.21) 

E[nl cg ) < 0.648 x 10" 2 ( 5 ) 2 

(7.22) 


The set of optimized feedback gains is shown in table 7.7. These gains differ signif- 
icantly from those shown in table 6.7 where the inner and outer loops are designed 
separately. Table 7.8 gives the damping ratio and natural frequency of the closed- 
loop system eigenvalues. A minimum damping of 0.65 is achieved through the use of 
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nonlinear constraints; this is an improvement over the minimum damping of 0.60 in 
the optimal outer-loop design in chapter 6. 

Figures 7.7 shows the command and broken-loop control frequency responses of 
the optimized design. Command and control-loop bandwidths for these designs are 
summarized in table 7.9. Note that the design has equal bandwidth in both the 
airspeed and altitude paths. The command bandwidth in the velocity variable is 
slightly less than the one of the outer-loop design in chapter 6, while the altitude 
command bandwidth is improved. Control-loop bandwidths are similar in the two 
designs. 

Figures 7.8 and 7.9 show time responses of the closed-loop systems to a step 
command in altitude of 1 000/2 and in airspeed of 10 fps respectively. Both command 
responses have almost the same time constant confirming the results of table 7.9. 

Table 7.10 summarizes the closed-loop rms responses to Dryden clear air turbu- 
lence at the 99% probability level of intensities. Results clearly illustrate the effective 
usage of direct bounds on the respective output covariances. By imposing inequality 
constraints on the control and output covariances, one is guaranteed at the outset 
that the optimal design, when converged, will have the same or lower control and 
output covariances than those of the classical design in table 5.4. 

Table 7.11 summarizes values of single- loop stability margins. Single-loop stability 
margins are found to be satisfactory and similar to those in table 6.1 1 of the previously 
optimized outer-loop design. 

Figures 7.10-7.12 show plots of minimum singular values of the return-difference 
and the inverse-return-difference transfer matrices respectively at both the control and 
sensor paths. In these plots, diagonal scaling on the transfer function matrices has 
been used to reduce conservatism and thereby improve estimates of actual multiloop 
stability margins. Table 7.12 gives guaranteed multivariable stability margins in 
both the control actuator and sensor paths according to equations (D.8) and (D.14) 
of Appendix D. Again these results are comparable to those in table 6.12 for the 
optimized outer-loop design of chapter 6. 

From the above results, one can deduce that to maximize tracking performance 
along with other requirements such as closed-loop damping, covariance responses to 
turbulence in an autopilot design problem, it is often necessary to re-examine the 
inner-loop design gains in conjunction with the outer-loop gains. We have demon- 
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strate in this chapter the usage of constrained nonlinear optimization in the synthesis 
of an airspeed and altitude-hold autopilot. Here the inner and outer-loop gains are 
optimized simultaneously to yield an overall improved design over those presented in 
chapters 5 and 6. 



Table 7.1: Optimal Inner- and Outer-Loop Feedback Gains (FLTl) 


Parameter 

SANDY 

I\tp 

0.56168 

Nep 

3.2326 

Kji 

0.33981 

Kei 

2.2854 

I< q 

3.2620 

Kb 

5.0000 

K v 

-0.010517 

K h 

-2.3882 xlO -4 

Kca is 

0.9322 

NfiW 

80,000 


Table 7.2: Closed-Loop System Poles (FLTl) 


Design 

SANDY 

Mode 

C 

io n (rad /sec) 

Phugoid mode 

0.8 

0.58 

Short period 

0.74 

3.13 

Altitude Mode 

1 

0.065 

Integral Speed Mode 

1 

0.125 


Table 7.3: Command and Broken-Loop Control Bandwidths (FLTl) 


Command/Control Path 

SANDY 


Bandwidths (rad/sec) 

h c 

0.06 

1 V C 

0.42 

Sec 

3.67 

Sthc 

0.40 





Table 7.4: Closed-Loop RMS Responses to Turbulence (FLTl) 


(<r u = 6.6 fps, a w = 6.3 fps) 


Variable 

SANDY 

7 (deg) 

0.9538 

V (fps) 

0.4755 

h (ft) 

13.1 

n *cg (g) 

0.07585 

Sth (lbst/lbsw) 

0.02718 

6' (deg) 

1.2982 


Table 7.5: Single- Loop Type Stability Margins (FLTl) 


Design 

SANDY 

Margins 

Gain Margin 
(dB) 

Phase Margin 
(deg) 

Actuator Paths 

(—35, -f-oo ) 

-57 

Sensor Paths 

(-19,-1-14.5) 

(-55, +70) 


Table 7.6: Guaranteed Multivariable Stability Margins (FLTl) 


Design 

SANDY 

Margins 

Gain Margin 
(dB) 

Phase Margin 
(deg) 

Actuator Paths 

(-14, +5.1) 
(-5.2, +14.5) 

±47.9 

Sensor Paths 

(-2.7, +3.95) 

±21 






Table 7.7: Optimal Inner- and Outer-Loop Feedback Gains (FLT2) 


Parameter 

SANDY 

K tp 

0.3965 

Kep 

6.549 

Kti 

0.27233 

Kei 

1.5469 

I< q 

5.3573 

K e 

-0.09335 

K v 

-3.593 xlO -3 

I<h 

-1.2038 xlO -4 

Kcas 

0.07006 

A GW 

80,000 


Table 7.8: Closed-Loop System Poles (FLT2) 


Design 

SANDY 

Mode 

C 

u> n (rad /sec) 

Phugoid mode 

0.65 

0.28 

Short period 

0.65 

3.27 

Altitude Mode 

1 

0.165 

Integral Speed Mode 

0.65 

0.16 


Table 7.9: Command and Broken-Loop Control Bandwidths (FLT2) 


Command/Control Path 

SANDY | 


Bandwidths (rad/sec) 

h c 

0.20 

v e 

0.20 

S ec 

4.0 

Sthc 

0.27 





Table 7.10: Closed-Loop RMS Responses to Turbulence (FLT2) 


(<7 U = 4.7 fps, <7u, = 4.7 fps) 


Variable 

SANDY 

7 (deg) 

0.1562 

V (fps) 

0.2059 

h (ft) 

5.5316 

n zcg (g) 

0.08050 

&th (lbst/lbsw) 

0.003878 

5 e (deg) 

0.1388 


Table 7.11: Single-Loop Type Stability Margins (FLT2) 


Design 

SANDY 

Margins 

Gain Margin 
(dB) 

Phase Margin 
(deg) 

Actuator Paths 

( — 17, -f-oo) 

-68 

Sensor Paths 

(-8.6, +12.7) 

(-50, +73) 


Table 7.12: Guaranteed Multivariable Stability Margins (FLT2) 


Design 

SANDY 

Margins 

Gain Margin 
(dB) 

Phase Margin 
(deg) 

Actuator Paths 

(-11, +4.7) 
(-5.7, +22) 

±55 

Sensor Paths 

(-2.9, +4.42) 

±23 
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Figure 7.1: Loop Frequency Responses of Optimal Inner and Outer-Loop Design 
(FLT1) 
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Figure 7.2: h c Command Responses of Optimal Inner and Outer- Loop Design (FLTl) 
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Figure 7.3: V c Command Responses of Optimal Inner and Outer- Loop Design (FLTl) 
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Figure 7.5: Robustness at the Control Paths of Optimal Inner and Outer- Loop Design 
(FLT1) 
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Figure 7.7: Loop Frequency Responses of Optimal Inner and Outer- Loop Design 
(FLT2) 
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Figure 7.8: h c Command Responses of Optimal Inner and Outer-Loop Design (FLT2) 
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Figure 7.9: V c Command Responses of Optimal Inner and Outer- Loop Design (FLT2) 
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Figure 7.10: Robustness at the Control Paths of Optimal Inner and Outer-Loop 
Design (FLT2) 



Figure 7.11: Robustness at the Control Paths of Optimal Inner and Outer- Loop 
Design (FLT2) 
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Figure 7.12: Robustness at the Sensor Paths of Optimal Inner and Outer-Loop Design 
(FLT2) 




Chapter 8 

CONCLUSIONS AND RECOMMENDATIONS 

8.1 Conclusions 

Traditionally, classical methods for control system design via successive single loop 
closure using root locus or Bode diagrams have been applied to the synthesis of 
feedback controllers for multiloop systems. Such designs often lack performance and 
robustness to uncertainties in the plant model. Automatic flight control systems 
(AFCS) on commercial aircraft are a prime example. Current AFCS design based 
on the single loop approach is evident in the architecture of the autothrottle and 
autopilot control systems. The autopilot is designed to regulate altitude through 
feedback to the elevator control surfaces while the autothrottle separately controls 
speed through feedback to the engines. 

The basic problem with classical design procedures is that they usually overlook 
the multiloop aspects of the system dynamics, neglecting significant cross-coupling 
among different feedback control systems. In the case of the AFCS, flight path and 
speed errors are produced by both the throttle and elevator controls such that closure 
of both the autothrottle and autopilot control loops may be destabilizing. Further- 
more, classical design procedures applied to multiloop controllers yield little insight 
into the problem of design tradeoff between performance and robustness to modeling 
uncertainties. Obtaining a satisfactory design with classical procedures is far from 
straightforward, and can be quite time consuming if previous design experiences are 
lacking. What is needed is a multivariable design procedure that directly incorporates 
design objectives in terms of closed-loop stability, performance, and robustness into 
the design tradeoff. 

Development in multivariable design techniques have concentrated on improving 
controller performance and robustness to modeling uncertainties. Modern control 
design methods such as LQG, H°°, and ^-synthesis allow designers to synthesize 
controllers that meet different measures of performance and robustness. Controller 
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designs obtained with modern techniques are usually of high order (i.e. the order 
of the controller is greater than or equal to the order of the plant model). Conse- 
quently, controller model reduction must be considered for practical implementation. 
Furthermore, design tradeoff between many real design goals (i.e. damping of system 
poles, disturbance rejection, command tracking, etc. ), while easier than with classical 
procedures, is still far from simple. 

In this study, the design method SANDY for robust controller design described in 
reference 6 is successfully applied to the design of an integrated autothrottle/autopilot 
control system. A complete design procedure encompassing key concerns such as 
closed-loop stability, control and command bandwidths, limited control activities to 
disturbance and command tracking performance has been developed for the synthesis 
of an integrated autothrottle/autopilot design based on the TECS concept. The 
usefulness of the procedure has been extensively demonstrated in actual design of an 
airspeed and altitude-hold system for a TSRV vehicle at two typical flight conditions. 
It should be noted that the constrained parameter optimization method implemented 
in the computer software SANDY has found wide application in aircraft flight controls 
and control of flexible mechanical systems, and provides a systematic approach for 
incorporating many real design objectives into the design tradeoff. Design goals 
such as damping of system modes, disturbance rejection, and command tracking 
are directly included into the design procedure. Direct tradeoffs are found to be 
systematic, and satisfactory designs are obtained with few iterations. 

The philosophy of total energy control (TECS) enables us to define a multiloop 
AFCS feedback control structure. Previous TECS designs use classical methods for 
obtaining the controller feedback gains and gain schedule. Results in this research 
demonstrate the application of the SANDY design algorithm to a practical flight 
control problem; design of an integrated autothrottle/autopilot control system. The 
procedure allows designers to address a number of objectives such as stability, distur- 
bance rejection and command tracking. Robustness to plant modeling uncertainties 
are maintained in the optimal design as compared with previous designs. 

Contributions of this research are in the following areas: 

• Formulation of performance indices for trade-offs in command and control band- 

widths. 
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• Innovative formulation of the nonlinear constraints for the closed-loop system 
eigenvalues. 

• Definition of a design variable V*(t) in terms of linear combinations of flight- 
path angle and velocity commands, and appropriate constraints on the design 
parameters to synthesize different levels of decoupling between the flight path 
command and velocity errors. This approach has a broad implication in the area 
of output decoupling control where cross-coupling between command paths can 
be systematically formulated. 

• Development of systematic design procedure for the TECS control laws with 
clearly defined design trade-offs using multiple performance indices, direct bound 
constraints, constraints on closed-loop stability and gain schedule across different 
flight condition. 

• Development of an integrated “inner loop” and “outer loop” design procedure 
for the TECS control-laws. 

8.2 Recommendations for Future Research 
Areas for future investigation include, 


• Examination of sensor implementation issues such as filtering of signals and 
estimation. 

• Extension of design procedures to directly include robustness measures into the 
design tradeoff. 

• Application of the design procedure to flexible-body aircraft with low structural 
damping. 

In this study, we have chosen to ignore the issue of sensor implementation for 
feedback of acceleration, flight path angle, and pitch damping variables. In practice, 
filtering of sensor signals is usually required. In fact, present TECS designs have 
included complementary filters for estimation of certain feedback quantities. The 
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design procedure developed in this investigation may be used to simultaneously design 
the feedback controller and the associated filter/estimator for the feedback sensors. 
Dynamic models of sensor signal conditioning can be incorporated into the synthesis 
model. The sensor filter/estimator parameters are then optimized along with gain 
parameters in the feedback controller. 

The design procedure outlined in this research allows for the direct tradeoff be- 
tween controller design goals of stability, performance, and robustness to modeling 
uncertainties. However, application of methods for incorporating robustness measures 
into the controller design procedure have not been extensively investigated. These 
techniques consist of including additional plant models for off-nominal conditions into 
the design procedure, and specification of H°° bounds on outputs of selected system 
transfer matrices for unstructured plant uncertainties. 

A useful extension of this research would be to investigate the design tradeoff with 
robustness to modeling uncertainties. Available techniques for analysis of structured 
singular values could be used in defining “worst-case” design conditions for inclusion 
as off-nominal plant models into the design procedure, as discussed in section 4.3. Fur- 
thermore, methods for improving multivariable stability margins through If 00 bounds 
on selected system transfer matrices could be explored. 

The design procedure developed in this research is well suited to the synthesis of 
controllers for flexible mechanical systems. An interesting application of this proce- 
dure would be to the design of an AFCS for a highly flexible aircraft, such as for 
modal suppression, gust alleviation, and ride quality control. Nonlinear inequality 
constraints for damping of system eigenvalues would be useful for ensuring minimal 
damping in the flexible-body modes. Numerical algorithms need to be extended to 
handle systems with repeated eigenvalues that are not diagonalizable. 
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Appendix A 


EVALUATION OF A QUADRATIC PERFORMANCE 
INDEX TO RANDOM INITIAL CONDITIONS AND 

IMPULSE INPUTS 

The finite-time performance index in equation (3.6) evaluated to random initial 
conditions or impulse inputs can be evaluated as in equation (3.7) for white-noise in- 
puts when certain conditions are met. This equivalence is established in the following 
sections. 

A.l Random Initial Conditions 

The closed-loop system associated with each plant condition has the following form 


x'(t) = A'x'(t) 

(A.l) 

y r (t) = Cpx'(t) 

(A.2) 

AO) = < 

(A. 3) 


where the initial conditions x'(0) are gaussian random variables with zero mean 
E\ x '<] = 0 a °d covariance E\x' 0 x'J] = X' 0 . The feedback control inputs can also 
be expressed in terms of the states of the closed-loop system as 

u p (t) = y e (t ) = C' c x'(t) (A. 4) 

The performance index evaluated to random initial conditions is 

M‘l) = \fj ^'E[y^(t)Qy r (t) + u T r (t)Ru r (t))dt (A. 5) 

Using equations (A. 2) and (A. 4), the performance index J\(tj) becomes 

M‘l) = {( c ? Qc; + Cf. RC' c ) jf 1 ' e , “ , £:[x'(()*' T (t)]<if } 


(A.C) 
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Solution of the closed-loop system responses to initial conditions is simply 

*'(<) = «' 1 V„ (A. 7} 

Combining equations (A. 6) and (A. 7) yields 

M*l) = jir {jf ' e< A ‘+°” Tl (c; T QC; + (7' r 72C')c*' 4 +0/ *'(fi A"} (A.8) 

The feedback controller must provide an asymptotically stable closed-loop system 
so that the limit in equation (A.8) exists and is finite. 

A. 2 White-Noise Inputs 

For random white-noise inputs, we consider the following closed-loop system 

x'{t) = A'x'(t) + r'w(t) (A. 9) 

Vp{ 0 = C' p x'{t) (A. 10) 

x'(0) = 0 (A. 11) 

where the disturbance inputs te(0 are gaussian random processes with zero mean 
£[te(0] = 0 and covariance £[te(<)u,’(r) T ] = W 0 6(t — r). The feedback control inputs 
can be written as 

MO = MO = c'x'(0 (A. 12) 

The finite-time performance index to random disturbance inputs is 

Mil) = (A. 13) 

where £J a ['] * s the expectation operator for a closed-loop system destabilized by a 
factor a. 

Substituting equations (A. 10) and (A. 12) into the expression for J 2 (</)> we obtain 
after some manipulations 

Mil) = \ir {(Cj r i QC' P + C? flC;)E„[ I '(( / )x' 7 '(i / ))} 


(A. 14) 
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where the time response x'(t) is given by 

x'{t)= f l e {A '+ aI * l - T) Y'w{T)dT (A. 15) 

Jo 

The performance index J 2 (tf) becomes 

Mh) = \ ir {{C?QC' v + C?RC' c ) (A. 16) 

f l} (*’ e {A ' +aI ^ t t-^r , E[w(T)w T (s)]r ,T e^ ,+aI)T ^- a UTds} 

Jo Jo 

Recalling that £'[u;(r)iu(.s) 7 ’] = W 0 8(t — 5 ), and with some rearranging J 2 (tj) becomes 

u * /) = c(A ‘ + °' ( ai7 > 

(C'J QC’ p + C?RC , c )e (A ' +oI){t i- T) dTT'W 0 Y n '} 

If we perform a change of variable t = tj — r, then 

Mh) = \Trace{j\( A '+°" Tt (A.18) 

(C?QC' p + Cf RC[)c (A '+ aI)t dtrW 0 T a '} 

Therefore it can be seen from equations (A. 8) and (A.18) that if 

T'W 0 T n ' = X' 0 (A. 19) 


then J\{tf) = J 2 (tj)- With this equivalence established, the problem associated with 
random initial conditions can be formulated as the one corresponding to random 
white-noise inputs. 

A. 3 Parameterized Random Impulse Inputs 

For random impulse inputs, the closed- loop system has the following form 


m 

= A'x'(t) + B'w(t) 

(A. 20) 

Vp(i ) 

= C' p x'(t ) 

(A.21) 

*'(0) 

= 0 

(A. 22) 
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where the inputs w(t) are parameterized random impulses i v(t) = ycmd$(t) with 
amplitude vector y^d being random with zero mean E[y cm d] = 0 and covariance 
E[y cm dyJmdl ~ Ycmd- The feedback control inputs can be written as 

MO = Vc(t) = C'x'it ) (A. 23) 

The finite-time performance index to random impulse inputs is 

J ^ tj ) = 2 L + Up( i )Ru p (t)]dt (A. 24) 

Substituting equations (A. 21) and (A. 23) into the expression for Jz(tj), we obtain 
after some manipulations 

Mtj) = ^tr^(C^QC; + C^RC:)j\ 2at E[x\t)x' T (t)}dt^ (A. 25) 

where the time response x f {t ) is given by the time-convolution integral as 

x'(t) = J‘ e ( A '+°m- r ) (A. 26) 

and since u>(r) = y cmd ^(r), we have 

i '(t) = J A '+°"‘B'y mi (A.27) 

The performance index Jz{tj) becomes 

Mh) = (A.2S) 

(C?QC; + C? RC'y A '*°"‘itE\B'y cmd y T md B' r \) 

Recalling that Ely^y^} = Y^, we have 

Mtj) = \tr | jT ' e (A ' +Ql Y\C^QC' p + C^RC , c )e (A ' +Ql)t dtB'Y cm dB rr yK.2d) 

It can be seen from equations (A. 18) and (A. 29) that if 

rW o r ,T = B'Y cmd B * (A. 30) 

then we have Jz(tj) = ^(t/). Thus, the problem associated with parameterized 
random impulse inputs can also be formulated as one corresponding to white-noise 
inputs if we use the equivalence relation of (A. 30). 


Appendix B 

CONTROLLER DESIGN FOR COMMAND TRACKING 
AND DISTURBANCE REJECTION 

Two common classes of control designs to external inputs are: command tracking 
or model-following and disturbance rejection. Control-law synthesis for these design 
problems can be formulated in terms of the minimization of a quadratic performance 
index of the form 

J = hm l ' E[yJ(t)Qy c (t) + uJ(t)Ru t (t)]dt (B.l) 

1 / JO 

The quadratic performance index J contains penalties on transient responses of 
the closed-loop system to random commands and/or disturbances. For command 
tracking, the variables y t (t) and u ( (t) in equation (B.l) are the tracking errors defined 

by 

y«(0 = y p (t)-y m (t) (B.2) 

«c(0 = Up(0 - Ucmd(t) (B.3) 

where y m (t ) are the outputs of a command generator model and u cmd (t) are the con- 
trol inputs that produce zero tracking errors in y c (t) at steady state. In order to 
ensure that the performance index in equation (B.l) be bounded, all the penalized 
variables in y e (f) and u ( (<) corresponding to the nonzero weighting matrices Q and R 
must converge asymptotically to zero as t — ► oo. For example, in the case of output 
commands y m (t) having nonzero constant steady-state values (i.e y m (t) = ycmd 0 n{t) 
where /z(<) is the unit-step), it can be shown through the set point calculations (sec- 
tion B.3) that the required controls u cm(i (<) at steady state are linear functions of 

Vcmda • 

For the problem of disturbance rejection, the tracking errors y t (t) and u t (t) in 
equation (B.l) are defined as 


y<(0 = y P {t) - y P „{t) 


(B.4) 
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u t {t) = u p (t)-u p „(t) (B.5) 

where y p „(t ) are the desired output responses while u Pat (t) are the necessary control 
inputs to maintain the desired outputs y p „{t ) in the presence of disturbances w(t). 
The relationships between u Paa (t), y Paa (t ) and the disturbances u>(<) are derived in 
section B.4. 

Command outputs y m (t) and disturbance outputs w(t) in the form of impulses, 
steps, ramps or a combination of these can be derived from linear time-invariant 
models with appropriately choosen initial conditions and impulse inputs. Details are 
given in sections B.l and B.2. 

B.l Models for Parameterized Random Commands 

Commands y m (t) in the form of impulses, steps or ramps can be generated using a 
linear time-invariant model of the form, 

im(t) = A m x m (t) + B m y cmd (t ) (B.6) 

y m (t) = C m x m {t) + D m y cmd {t) (B.7) 

with x m (0) = x mo and y cmd (t) = y cm do^(i)- The quantities x mo , y cmdo along with the 
state matrices A m , B m , C m and D m define the desired outputs y m (t). 

For example, a command vector of impulses can be created using 

1/m ( 0 l/cmd o ^(0 (B.S) 

where the state vector x m (t ) is of zero dimension and the matrix D m = I. The 
impulse magnitude ycm do is a vector of gaussian random parameters with zero mean 
E[ycmdo] = o and covariance = Fcmdo • 

On the other hand, a vector of step commands can be obtained from 


*m(<) 

= (0)x m (<) + ycmdoHi) 

(B.9) 

S/m(0 

= Xm{t) 

(B.10) 

x m (0) 

= 0 

(B.l 1) 
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yielding a step command vector y m (t) = where y.(t) is a unit-step function 

an d y C md 0 is a vector of gaussian random parameters with zero mean E\y cmdo \ = 0 
and covariance E[y cmdo yJ mdo ] = Y^. 

An alternate formulation for the vector of step commands is 


im(t) = 

(0)* m (f) 

(B.12) 

y m (t) = 


(B.l 3) 

*m(0) = 

X mo 

(B.l 4) 


In this case, the resulting step command is y m {t) = x mo n(t) where x mo is a vector of 
gaussian random parameters with zero mean £[!„,] = 0 and covariance £’[i mo x^ ] = 
* mo . 

Filtering of command outputs can be done by introducing additional dynamics 
into the command generator model described in equations (B.6)-(B.7). The choice of 
the state matrices A m , B m ,C m and D m will completely define the command signals 
J/m(0- For example, if the system matrix A m has a pair of eigenvalues at A = ±ju \ 0 
(i.e on the imaginary axis), then the outputs y m (t) will have components of sinusoidal 
functions of frequency u> 0 . Figures (B.l) and (B.2) show time responses of a stable 
first-order filter to a scalar impulse command and a step command respectively. 

B.2 Models for Parameterized Random Disturbances 

In a similar manner, external random disturbances w(t) in the form impulses, steps 
and ramps can be generated using a linear time- invariant model of the form, 

id(t) = A d x d (t) + B d w d {t) (B.15) 

w(i) = C d x d (t) + D d w d (t) (B.16) 

with i<f(0) = x d(j and w d (t) = w 0 6(t). The quantities w 0 along with the state 
matrices A d , B d , C d and D d are parameters that define the disturbances w(t). 

For example, a disturbance vector of impulses can be created using 


iw(t) = w 0 8(t) 


(B.17) 
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where the state vector x d {t) is of zero dimension and the matrix D d = I. The 
amplitude w 0 is a vector of gaussian random parameters with zero mean £[u> 0 ] = 0 
and covariance E[w 0 Wq ] = W 0 . 

On the other hand, a vector of constant disturbances can be obtained from 


id(t) 

= (0)id(f) -f w 0 S(t) 

(B.1S) 

w{t) 

= x d (t) 

(B.19) 

Xd{ o) 

= 0 

(B.20) 


yielding a disturbance vector w(t) = w 0 p(t) where w 0 is a vector of gaussian random 
parameters with zero mean .Efieo] = 0 and covariance E\wqWq ] = Wq. 

Note that an alternate formulation for the vector of constant disturbances can be 
done using random initial conditions as described below, 

id{t) = (O)xj(O (B.21) 

w(t) = x d (t) (B.22) 

3M(0) = Ido (B.23) 

In this case, the resulting step disturbance vector is u>(<) = i d Q p{t) where x do is a 
vector of gaussian random parameters with zero mean E[x do \ = 0 and covariance 

£[Vd 0 J = to- 
other parameterized random disturbance functions can be produced, similarly to 
the command generator models of section B.l, by introducing additional linear filters 
into the disturbance model of equations (B.15)-(B.16). 

Random stochastic processes can also be modeled using equations (B.15)-(B.16) 
where, in this case, the disturbance inputs Wd{t) are white-noises with zero mean 
E[w d {t)} = 0 and covariance E[w d (t)wJ (t)] = — r). Outputs w(t) of the dis- 

turbance model are depicted as colored noises having power spectral characteristics 
defined by the disturbance state matrices A d , B d , C d and D d . 

B.3 Steady-State Responses to Command Inputs 

Command tracking and model-following designs involve the formulation of a perfor- 
mance index that contains tracking errors y c (t) = y p (t) —y m (t) and control feedback 
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errors u ( (t) = u p (t) — Ucm^t) as depicted in equation (B.l). These error variables 
can always be written as outputs of a closed-loop system. The quantities y m (t) are 
commanded outputs while u cm d(t) are the control inputs used to maintain zero steady- 
state tracking errors in y c (t). We would like to establish the necessary relationships 
between the commanded outputs y m (t) and the control inputs Ucm<i(0- 

Let’s consider a linear time-invariant plant model with control inputs u p (t) and 
outputs y p (t), 


x p (t) = A p x p {t) + B p u p (t) (B.24) 

y P (t) = C p x p (t ) + D p u p {t) (B.25) 


In the following derivation, we restrict ourselves to the case where the commanded 
outputs y m (t ) have only constant steady-state components; that is, we exclude any 
commanded outputs that are sinusoidal or grow unbounded with time such as ramp 
functions. This assumption enables us to treat this problem as one corresponding to 
the calculation of setpoints. Steady-state system responses are determined by letting 
i p (0 = 0, u p (t) -* u cmd (t) and y p (t) -> y m (t) in the limit as t — ► oo. Dimension of 
y P (t) (and thus y m (t)) must be equal to the dimension of u p (t) so that we have 


or 



A p B p / 0 \ 

C V D P \ y p {t ) ) 


(B.26) 


( * p ..(0 ) 

\ u cmd{t) J 


A p B p I 0 \ 

C P D p \ y m (t) ) 


(B.27) 


where the notation (•)„ denotes the quasi-steady values. We should emphasize 
that these relationships are correct only when the steady-state conditions have been 
reached. Nonetheless one can use relations in equation (B.27) to express the desired 
states and control command inputs as a function of the commanded outputs y m {t) 
for all time t > 0. The control command inputs «„„<<(*) and the state vector x p „(t) 
are therefore functions of time since the commands y m (t) may be time-varying, and 
they will settle to constant values as t — * oo. Note that the system matrix formed 
by the quadruple A p , B p , C p and D p in equation (B.27) is square and it must also 
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be invertible for the setpoint calculation. The quantities y m (t ) are outputs of a com- 
mand generator model and the values x Pst (t) and Ucmd(0 are computed according to 
equation (B.27). 

Let’s denote the inverse of the system matrix as 

A p B p 

_ C p D p 

where the matrix partitions K through N have the same dimensions as their coun- 
terparts A p through D p on the left-side of the equation. With these definitions, the 
command inputs Ucm^t) in equation (B.27) can be written explicitly in terms of y m (t ) 
as u cmd(0 = A f j/ m (<). The plant synthesis model in equation (B.24) for a command 
tracking problem is given by, 


i p (0 = 

A p x p (t) + B p u p (t ) + (0 )y m (t) 

(B.29) 

z/<(0 = 

CpXp(t) + D p u p (t) - Iy m (t) 

(B.30) 

u c (t) = 

(0)xp(0 + Iu p (t) - Ny m (t) 

(B.31) 


K M 
L N 


(B.28) 


where y m {t) is the command input vector discussed in section B.l. The tracking 
errors y e (<) and u ( (t) in the performance index of equation (B.l) are now outputs of 
the plant model in equations (B.29)-(B.31). 

To complete the problem formulation, plant model in equations (B.29)-(B.31 ) must 
be augmented with a command generation model as described in section B.l to form 
the final synthesis model in a command tracking problem. 

B.4 Steady-State Responses to External Disturbances 

The design problem for disturbance rejection involves the formulation of a perfor- 
mance index as in equation (B.l) that contains tracking errors y { (t) = y P (t) — y p ,,(t) 
and control feedback errors u e (f) = u p (t) —u Ptt (t). These error variables can be writ- 
ten as outputs of a closed-loop system. The quantities y p „(t) are the desired outputs 
of the closed-loop system in the presence of disturbances. The control inputs u p „(f) 
are used to reject the unwanted disturbances in the output responses and thereby 
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enforce zero tracking errors in y c (t) at steady state. We need to establish a set of re- 
lations similar to equation (B.27) of the tracking problem for the desired state vector 
x p „(t) an d control inputs u Pa$ (t). 

Let’s consider a linear time-invariant plant model with control inputs u p (<), dis- 
turbance inputs w(t) and outputs y p (i), 


•*-p(0 — A p x p (t) + B p u p (t ) + I>(<) (B.32) 

y p (t) = C p x p (t ) + D p u p (t) + Q p w(t) ( B.33 ) 


As before, the desired steady-state system responses are determined by letting x p (t) = 
y P (t) — ♦ y P „(0i u p(0 — 1 ► Wp,.(0 and w(t) — > to„(<) in the limit as t — ► oo. Again, 
the dimension of y p (t ) must be equal to the dimension of u p (t) so that we have 


/ x p (t) 
\ u p(0 



Ap B p I -T p w{t) \ 

C p D p _ \ yp(t) ~ ttp w (t) ) 


(B.34) 


where the notation (•)„ again denotes the quasi-steady values. The control input 
u p».(0 a °d the state vector x P)$ (t ) are now functions of the external disturbances 
anc ^ the desired outputs y p „{t). And they may be time-varying. 

Note that the system matrix formed by A p , B p , C p and D p in (B.34) is a square 
matrix and it must have full rank for its inverse to exist. 

We denote the inverse of the system matrix as 


(B.35) 


' A* 

B P ' 

-1 

’ I< 

M ' 

. c. 

Dp. 


L 

N 


where the matrix partitions K through N have the same dimensions as their counter- 
parts A p through D p on the left-side of equation (B.35). The plant synthesis model 
in equation (B.32) for a disturbance rejection problem is given by 


x p {t) = A p x p (t) + B p u p (t) + (0)y p Jt) + r p w„(t) (B.36) 

y c (t) = C p x p (t) + D p u p (t)~ Iy p Jt) + n pW , 3 (t) (B.37) 

“«(0 = (Q)xp(t) + Iup(t)-Ny p „(t) + Nn p w S 3 (t) (B.38) 


The tracking errors y ( (t ) and u { (<) in the performance index given in equation (B.l) 
are now outputs of the plant synthesis model in equations (B.36)-(B.38). 
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To complete the problem formulation, the plant model in equations (B.36)-(B.38) 
must be augmented with a disturbance model as described in section B.2 to form the 
final synthesis model in a disturbance rejection problem. 



Figure B.l: Response of a First-Order Stable Filter to an Impulse Command 


Commanded Output 
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Appendix C 


NONLINEAR CONSTRAINTS ON SYSTEM 

EIGENVALUES 


Design requirements involving damping ratio and frequency of closed-loop eigen- 
values are one of many important design criteria that a feedback control system must 
satisfy. For example, flying qualities in the longitudinal axis are often expressed in 
terms of desired locations of the short-period and phugoid modes. To achieve the 
desired characteristics in the closed-loop eigenvalues, one common design procedure 
based on optimal control would be through the use of output penalty variables in the 
quadratic performance index. Clearly, a design procedure for closed-loop stability via 
a performance index is indirect and necessitates numerous trials before arriving at a 
reasonable set of penalty variables. 

The reason behind this difficulty is that there is no one-to-one correspondance 
between location of closed-loop system eigenvalues and transient responses of output 
variables to disturbances. A more direct approach for the assignment of closed-loop 
eigenvalues is with the use of direct constraints. The following sections will describe 
the formulation of two types of eigenvalue constraints: one type of constraint is on 
the real parts of the eigenvalues, and the other is on their damping ratio parameters. 
Other variation or extension from these set-ups can be easily implemented by a simple 
modification of the proposed schemes; for example, the case involving damping and 
frequency constraint requirements that are frequency-dependent. The advantage of 
the proposed formulation of eigenvalue constraints is that the constraint functions 
in both cases are continuously differentiable with respect to the design parameters 
in the system matrix. However the above statement holds only for systems that are 
diagonalizable, i.e having a non-degenerate Jordan structure. 


C.l Constraints on the System Eigenvalues 

Consider a linear time-invariant system described by the following state model 
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i(t) = Ax(t) (C.l) 

where x(t) is a state vector of dimension n, and A is a system matrix of dimension 
n x n. System stability is governed by the eigenvalues of the system matrix A defined 
by 

Avi = A,u,, (i = 1, n) (C.2) 

where A, is the i th eigenvalue of the system matrix A and v, is the corresponding 
eigenvector. For practical purposes and from here on, we assume that the system 
matrix A is diagonalizable, i.e. there exists a nonsingular similarity transformation T 
constructed from a set of n completely independent eigenvectors [w, •],•_, n such that 

T-'AT = A (C.3) 

= diag (A,) t=1 „ (C.4) 

In this appendix, we examine constraint formulation on two basic parameters of an 
eigenvalue A, = <r, + (i = l,...,n) where j = y/C- l 5 namely: (1) its real part cr, 
and, (2) its damping ratio (, defined as 


C. = 



(C.5) 


Clearly for stability, one must have at least o x < 0 or (, > 0, (i = l,...,n). Usually 
requirements for closed-loop stability are expressible in terms of desired real parts 
and damping of system eigenvalues as indicated in figure C.l. 

The corresponding constraint inequalities are 


^ (* — 1 , ..., 71 ) (C.6) 

C. ^ Cmin) = 1 ,..., n) (C.7) 

where <r m ax a nd are design specification parameters. Equations (C.6) and (C.7) 
constitute a set of 2n inequality constraints on the system eigenvalues. For conve- 
nience and simplicity of implementation, one can employ an equivalent formulation 
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that yields basically the same results, but would involve only a single inequality con- 
straint. To understand the achieved accuracy level of a constrained eigenvalue under 
a given tolerance on the inequality constraint, we need to examine, as given in the 
following sections, the results corresponding to the constraints on its real part and 
its damping ratio separately. 

C.1.1 Constraints on the Real Part of System Eigenvalues 
An equivalent formulation of equation (C.6) is as 

0 < <7 < \-t 2 (C.S) 

where c is a preselected tolerance level and the variable a is defined as 

1 71 

° = o J2 [ ma *(<A - <7mi»x,0)] 2 (C.9) 

1 . = 1 

with a m a X being the desired maximum real part of the eigenvalues. Note that eigen- 
values that satisfy the constraints in equation (C.6) do not contribute to the variable 
tr; only those eigenvalues that violate the constraints are accounted for in equation 
(C.9). Clearly, when the constraint in equation (C.8) is satisfied then the real parts 
of the active eigenvalues cr,(i = l,...,n) (i.e. those that contribute to a in equation 
(C.9)) are within a tolerance ±e from the desired bound < 7 max . Note that the param- 
eter a is a continuously differentiable function with respect to the design parameters 
in the system matrix A. Section C.2 derives the sensitivity function with respect 
to a parameter p in the system matrix A that is diagonalizable. 

C.1.2 Constraints on the Damping Ratio of System Eigenvalues 
An equivalent form of equation (C.7) is 

0<C<^ 2 (C.10) 

where e is a preselected tolerance level and the variable £ is defined as 

1 n 

C = ^E [max(<7,+ | w, | cota,0)] 2 

£ .=i 


(C.ll) 
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with cos cr = £„„„ and sin a = y 1 — where is the desired minimum damping 
ratio applied to all the system eigenvalues. Physically the constraint variable £ is the 
sum of the square of the horizontal distance d, for the eigenvalues with nonnegative 
imaginary parts (i.e. uj, > 0) that are located above the damping line (figure C.2). 
Any eigenvalue with nonnegative imaginary part that falls below the damping line 
does not contribute to the constraint variable £. Those eigenvalues with negative 
imaginary parts (i.e. u){ < 0) will contribute the same amount as their complex 
conjugate counterparts. When the constraint in equation (C.10) is satisfied, the 
damping ratio £, of each active eigenvalue A ,( i = l,...,n) (i.e. one which contributes 
to £ in equation (C.ll)) will be bounded by 



Figure C.3 shows relative levels of accuracies (i.e values of £/ £ mm — 1 ) achieved from 
the constraint equation (C.10) as a function of the ratio j^-r, i.e. plots of 8± £min) • 
Thus for an actively constrained eigenvalue A,, the desired damping ratio £,„;„ is 
satisfied to within an acceptable level of accuracy, say 0.01, if the parameter t is 
selected such that 

6 

To ensure that all eigenvalues A,(i = l,...,n) have the desired damping ratio £m] n 
to within a relative accuracy of at least 0.01, the parameter e must be chosen based 
on the minimum value of | cr, | (i = l,...,n). Together with the constraint equation 
(C.8) for the real part of the eigenvalues, a reasonable estimate for the minimum of 
| <t, | (t = 1, ...,n) would be | cr max | . 


°i I 


,c- 


<0.01 


(C.15) 
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C.2 Gradients of the Constraints with Respect to Parameters in the System Matrix 

To apply the constraints of equations (C.8) and (C.10) into a control design prob- 
lem using parameter optimization, it is often required that gradients of the con- 
straint functions (e.g. o and £) with respect to a system matrix parameter (e.g. 

element (fltj)ij-, „) be supplied to the nonlinear optimization algorithm (in this 

case, NPSOL [Ref.8]). Determination of gradients using finite- difference approxima- 
tion can be time-consumming and numerically inaccurate. Whenever possible, it is 
more efficient and accurate if expressions of constraint gradients can be derived and 
evaluated analytically. 

Recall that constraint on real part of the system eigenvalues is formulated as a 
constraint on o (equation (C.9)) where 

1 " 

a = nE[ max ( ff '- -^n»x,°)] 2 (C.16) 

1 .=1 

Gradient of a with respect to a parameter p in the system matrix is given by 

da " d*i 

Qp / . max((7, ^inan 0 ) (C.l i) 

Constraint on damping of the system eigenvalues are formulated as a single constraint 
on ( (equation (C.ll)) where 


(= nE [max(fT I + | u/ t - | cot o, 0)] 2 


i=i 


(C.1S) 


Gradient of £ with respect to a parameter p in the system matrix is given by 

n ( da 8 \ 

- E max (*.+ I I cot a, 0) + S g n ( W| . ) co t a J (C.19) 

where sgn(-) is the sign function (i.e. sgn(x) = 1 if x > 0, sgn(x) = 0 if x = 0, and 
sgn(x) = -1 if x < 0). 


C.2.1 Sensitivity Function of System Eigenvalues 

In this section we derive analytically the sensitivity function of the system eigenvalues 
A,, (f = l,...,n) with respect to a parameter p in the system matrix, i.e. 
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From the definition of system eigenvalues given in equation (C.2) we have 


Av, = A (i = 1, n) 


(C.20) 


Differentiating both sides of equation (C.20) with respect to a parameter p (where p 
is an element n of the system matrix A) to yield 


dA dvi d \ , dv, 

~E~ v i + -4—— = — — Vi + A,' — 
op op op op 


( 0 . 21 ) 


or 


/, \ T \ dv i dX ' dA 

(/l - k,) T P ~ = -fy Vi 

Let T be a transformation that diagonalizes the system matrix A as 


(C.22) 


T~ l AT = A 


= diag{\i) i=l n 


(C.23) 

(C.24) 


Multiplying both sides of equation (C.22) by T 1 and after some manipulations we 
have 


A x 0 • • • • 0 

0 • 

■ 0 A,_i 0 

0 0 0 

0 A 1+1 0 

• 0 

o • • • • 0 A n 


^-x dv, , d \ , 

T st- T vi w 


T - A 

= -T -t-v, 

dp 


(C.25) 


Clearly, one can solve for the sensitivity ^ from the i tfl row of equation (C.25) as 


ax, (r-'tv,),. 

dp (T-'»,), 


(C.26) 
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where the notation (•); denotes the i th row of the enclosed quantity. Recall that 
the i lh eigenvalue is A, = a, + ju t , thus the sensitivities ^ and ^ are simply 
the real and imaginary parts of ^ respectively. Note that the computation of the 
sensitivity function involves the determination of a transformation T consisting of 
system eigenvectors v,(i = For the case where some of the eigenvalues are 

repeated and the system does not have a complete set of eigenvectors corresponding 
to these eigenvalues, then the sensitivity function associated with these eigenvalues 
does not exist. 

For a feedback controlled system, one is usually interested in obtaining the sen- 
sitivity of the closed-loop eigenvalues with respect to the controller design gain pa- 
rameters. In the next section we will establish the relation between the closed-loop 
system matrix and the controller state matrices. 

C.3 Sensitivity of the Closed-Loop System Matrix with Respect to Controller State 

Matrices 

Consider a plant model with the following state equations 

i P (l) = A p x p (t) -f B p u p (t) (C.27) 

y p (t) = CpXp(t) + DpUp(t) (C.2S) 

where x p {t) is the plant state vector of dimension n, u c (t) the control input vector 
of dimension m and y p (t) the sensor output vector of dimension p. A linear time- 
invariant feedback controller for the above plant is of order r (where r < n) and has 
the following state-space description, 

*e(0 = A c x c (t) + B c u c (t) (C.29) 

Vc{t) = C c x c (t) + D c u c (t) (C.30) 

where x c (t ) is a controller state vector of dimension r and (A C ,B C , C c , D c ) are the 
controller state matrices of appropriate dimensions. Note that, in a feedback config- 
uration, the plant sensor output vector y p (t) is connected to the input vector u c (t) 
of the controller, and the outputs y c (t) of the controller are the inputs u p (i) of the 
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plant. For convenience, we define a matrix K that contains in a compact form the 
quadruple (A c , B c ,C c , D c ) as follows, 


D c C c 
B c A c 


(C.31) 


In this formulation, when we refer to an element of the matrix K, we actually address 
one of the elements of the controller state matrices A c , B c ,C c , and D c . 

With the feedback control loops closed and assuming that D c D p = 0, we obtain 
the following state model for the closed-loop system, 


/ x p (t) 

\ i c (0 


A p + B p D c Cp BpC c 
B c (I + D v D c )Cp A c + B c D p C c 



x p (t) \ 

X c {t) ) 


(C.32) 

(C.33) 


where A c i 0 , e d is the closed-loop system matrix. Stability of the closed-loop system 
is determined by the eigenvalues of the matrix A c i oae d. The sensitivity of the matrix 
Adoacd with respect to a parameter in the controller state matrix K (equation (C.31)) 
can be easily computed if we rewrite the matrix A c i oae d as follows, 


Acioaed — Ao A ( Bq + To A Do) A Co (C.34) 

where the matrices Ao, B 0 , Co, D 0 and 7 0 are independent of the controller matrix K 
and they are given by 


A o 
Bo 
Co 
D 0 


A P 

0 

B v 

0 

0 

D P 

0 


0 

0 

0 

I 

0 

I 

0 

0 


J (n+r)x(n+r) 


J (n + r) x (m+r) 


J (p+r)x(n+r) 


J (p+r)x(m + r) 


(C.35) 

(C.36) 

(C.37) 


(C.38) 
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h 


0 0 
0 I 


(n+r)x(m + r) 


(C.39) 


Let p be the ij th element of the controller matrix K (i.e. p = A*, ) then the sen- 
sitivity dA c i 0 , ed /dp is obtained by differentiating equation (C.34) with respect to p 
as 


dA 


closed. 


dp 


= {B 0 + r 0 A tJ D 0 ) KC 0 + {B 0 + I 0 KD 0 ) A U C 0 


(C.40) 


where A tJ is a matrix of dimension (m + r) x (p + r) with zero entries except for the 
ij element where it is equal to unity. Finally sensitivity of the closed-loop eigenvalue 
with respect to a controller parameter p can be obtained using equation (C.40) and 
equation (C.26) for the system matrix A c i ose j. 

Results presented in this appendix have been implemented in the design algo- 
rithm SANDY which combines the two types of nonlinear constraints given in equa- 
tions (C.8) and (C.10) into a single inequality constraint. Precisely, the eigenvalue 
constraint is defined as 


1 9 

-oo < c A < -e 2 (C.41) 

where c A is given by c x = \ £ df with 

»=1 

f max(<7,+ | u, I cota,0) ,x,>0 

a, = < ((j 

{ max{<j x -<r mAX ) ,x t - < 0 * 

and x, — cr max + j u;, | cota. Design specifications on the system eigenvalues are 
defined in the parameters <r max > 0 and 0 < a < t/2 where cosa = ( m ,' n . 
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Figure C.l: Constraint Boundaries for System Eigenvalues 



Figure C.2: Definition of Eigenvalue Real-Part and Damping Ratio Constraints 






Appendix D 

MULTIVARIABLE ROBUSTNESS MEASURES 


Multivariable stability margins are defined in terms of maximum allowable varia- 
tions of gain and phase in the feedback loops according to figure D.l for multiplicative 
type uncertainties. The complex gain matrix K is given by 


I< = (k 1T 




) 




(D.l) 


where m is the number of feedback controls. Note that K is not a function of the 
complex frequency s. The border-line for closed-loop stability is when the return- 
difference transfer matrix of a nominally stable closed-loop system becomes singular 
at some frequency 3 = jco, i.e. 


det [I + G(jtj)K) = 0 


(D.2) 


DA Robustness Test I 

For multiplicative type uncertainties as depicted in figure D.2, instability occurs when 
the return-difference matrix becomes singular for some value of s = juj y i.e. 

det [I + G(ju){I + L(ju>))] = 0 (D.3) 

or 

det = 0 (D.l) 

Closed-loop stability is therefore guaranteed when 

*Mw))<2.{l + G- l (ju>)) (D.5) 

for u > 0. The operators a(-) and ct( ) denote the minimum and maximum singular 
values respectively. In terms of the complex gain matrix I< of equation (D.l), the 
above bound becomes 


a(I< - I) < a (l + G~ l (ju)) 


(D.6) 
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If the complex gain matrix is restricted to a diagonal matrix of the form 

I< = diag (D.7) 

V / 7=1,..., 771 v ' 

then the inequality bound reduces to 

™< x m \/(1-M 2 + 2A;, 9 (1 -cos0„) < mma(l + G~'(ju;)) (D.8) 

With this relation, gain and phase margins for type I uncertainties are easily computed 
from a plot of minimum singular values of the inverse return-difference matrix of the 
nominal system. Note that this test requires that the nominal loop gain matrix G(jtj) 

be nonsingular. If G(ju>) becomes singular at some frequency w then this robustness 
test is not possible. 

D.2 Robustness Test II 

Tor multiplicative type uncertainties as depicted in figure D.3, instability occurs when 
the return-difference matrix becomes singular for some value of s = ju>, i.e. 

det [I + G(ju>)L(ju>)] = 0 (D.9) 

or 

det [/ + G(jw) + L~ l (ju) — /J = 0 (D.10) 

Closed-loop stability is therefore guaranteed when 

a(L-'(ju)-l) <<r(I + G(ju)) (D.l 1) 

for u) > 0. In terms of the complex gain matrix I< of equation (D.l), the above bound 
becomes 

_/) < 2.(1 + G(ju>)) (D.12) 

If the complex gain matrix is restricted to a diagonal matrix of the form 

K = diag m 


(D.13) 
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then the inequality bound reduces to 

!™ x m \ f l ~^ )2 + ^; (1 - cos ^) ^ m j n ^( 7 + (d 14 ) 

With this relation, the gain and phase margins for type II uncertainties are easily 
computed from a plot of minimum singular values of the return-difference matrix of 
the nominal system. 

Both equations (D.8) and (D.14) will yield bounds which are applicable to simul- 
tanous gain and phase variations in every control loop. Stability margins associated 
with each test may be combined to produce less conservative bounds. For example, if 
the smallest a (I + G(ju>)) is 0.7 and a (/ + G -1 (ju>)) is 0.6, then in the presence of 
phase uncertainties of ±30 degrees, the closed-loop system will tolerate simultaneous 
gain changes of —2.64 dB to 8.49 dB and —5.44 dB to 1.57 dB respectively. Remem- 
ber that robustness test I requires a nonsingular loop system matrix G{ju>). If this 
is not the case, then robustness test I can not be used to improve the estimates of 
multivariable stability margins. 

D.3 Effects of Scaling on Robustness Bounds 

Estimates of multivariable stability margins of type I and type II computed from 
equations (D.8) and (D.14) can be improved by scaling the system matrices G{ju), 
since singular values are not invariant under a similarity transformation. However, 
the characteristic equation of the nominal closed-loop system is unaffected by scale 
changes. This follows from the fact that, while det (/ + RG(jcj)R“*) = 0 is equivalent 
to det (/ + G(jij)) = 0 for any nonsingular transformation /?, their singular value 
norms <r (/ + RG(ju>)R-') and a (I + G(ju>)) can be much different. 

In fact, bounds on gain and phase variations determined from equations (D.S) 
and (D.14) can be improved by rescaling the loop transfer matrix G(s). If the type of 
loop uncertainty is restricted to a diagonal complex gain matrix in equation (D.7) and 
we only consider diagonal scaling matrices D — diag(d,) l=lt m , then the left-hand 
side of equations (D.8) and (D.14) remain unchanged from the scaling of the return- 
difference and inverse return-difference transfer matrices. Increasing the minimum 
singular values through diagonal scaling will give less conservative estimates of the 
multivariable stability margins determined from equations (D.S) and (D.14). 
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DA Osborne’s Method for Pre-Conditioning of Matrices 

Reference 22 describes a method for pre-conditioning of matrices. The procedure 
involves the minimization of the euclidean (or Frobenius) norm of a square matrix 
through diagonal similarity transformations. The method is easily programmed for 
numerical solutions and is seen to converge quite rapidly. The euclidean norm of the 
matrix A is defined as 

M’ = EEKf ( D15 > 


i=i j=i 


The minimization problem becomes 


min \\DAD~' ||* 


(D.1C) 


over all real diagonal scaling matrices D. For simplicity, we can take each diagonal 
element d } of D = diag(di, ...,d n ) as a free variable. The product of DAD~ l is 




(0.17) 


Thus, as the element dj varies, only the j th column and j th row of the matrix DAD~ A 
are affected, leaving the jth diagonal element of the product unchanged. The square 
of the euclidean norm \\DAD~ l ||^ can be rewritten to explicitly show the dependency 
on the single varying parameter dj as 


n ” (Pa 2 
r' < u 


I dad-x = 

t=l i=l a J 


(D.18) 


= + EE^ + MMD.19) 

J= . a ] a k ! I .,1 I 

}?k \ijtk j^k / 


= d 2 k R+—S + C 


(D.20) 


This expression is minimized when 


(D.21) 
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Therefore the euclidean norm of A may be minimized with respect to the element 
dk of D in one step. Osborne’s algorithm consists of stepping through each element 
d k from k = and repeating the process over again until the norm [jZ?y4Z) — 1 || e 

can no longer be improved and therefore has converged to the global minimum. It 
should be noted that the method is guaranteed to work only if the matrix A is in 
a form that is irreducible by permutation matrices. Otherwise, the algorithm would 
fail with either quantities S or R being equal to zero. 

D.5 Improving the Multiloop Stability Margins through Diagonal Scaling 

Our problem for improving the estimates of multivariable stability margins with di- 
agonal scaling is stated as 


max£ (DU + G(jw))D-') 

Recall that the minimum and maximum singular values are related by 

= i(^y 

so that our problem can be reformulated as 


(D.22) 


(D.23) 


min a- (£>(/ + G(ju>))~ l D~ l J (D.24) 

The maximum singular value of a matrix A is related to the euclidean norm of A 
as 

< Mile (D.25) 

Hence, minimizing the euclidean norm of D(I + would lower the upper 

bound on d ( D(I + G{ju))~ x D~ l ) in equation (D.24). So if the loop transfer matrix 
is poorly scaled, Osborne’s method provides a way for improving the estimates of the 
actual multivariable stability margins. 
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Figure D.l: Feedback System for Evaluation of Multivariable Stability Margins 



Figure D.2: Type I Multiplicative Perturbation where L(s) = 0 (nominal) 



Figure D.3: Type II Multiplicative Perturbation where L(s) = I (nominal) 




Appendix E 

TSRV B-737 DYNAMIC MODELS 

Longitudinal aircraft rigid-body dynamic models are obtained for the NASA Lan- 
gley Transport Systems Research Vehicle (TSRV). It is a highly modified Boeing 
737-100 aircraft for research in advanced guidance, navigation, control and display 
concepts [Refs. 12-14, 24]. Linear dynamic models in the wind-axis are given in the 
form 


— ApXp{i) + BpUp(t) 

(E.l) 

State and control input vectors are 


I p(0 = {«(0.“(0.«(0.*(0} 

(E.2) 

and 


<« = {*r(0, *.(<)} 

(E.3) 

respectively where 



u(t) = forward velocity component (fps) 

ci(t) = angle of attack (rad) 

q(t) = pitch attitude rate (rad/sec) 

0(t) = pitch attitude (rad) 

S T = total engine thrust (lbs) 

S e — elevator control position (deg) 

Two flight conditions were chosen for this investigation. One is an approach to landing 
condition and the other is a high-altitude cruise flight condition. Parameters for each 
flight condition are listed in tables E.l and E.2 for trimmed flight. The state dynamic 
matrix A p and control input matrix B v for flight condition FLTl are 

-0.4697 x lO" 1 17.09 0 -32.12 

-0.1524 x 10 -2 -0.7136 0.9995 0.8295 x 10~ 2 

-0.4397 x 10" 3 -1.235 -0.5020 -0.3233 x 10" 3 

0 0 1.0 0 
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0.40177 x 10 -3 

0 

-0.89059 x lO" 7 

-0.77001 x lO -3 

0.63002 x lO" 5 

-0.19737 x lO" 1 

0 

0 


(E.5) 


For flight condition FLT2, the matrices A p and B p are 


A p 


B p 


-0.7055 x 10~ 2 
-0.1275 x 10 -3 
0.3532 x 10" 4 
0 

0.40210 x 10~ 3 
-0.10185 x lO" 7 
0.62982 x 10" 5 
0 


21.30 0 

-1.032 0.9976 

-6.948 -1.009 

0 1.0 

0 

-0.10562 x lO" 2 
-0.99133 x 10" 1 
0 


32.17 

0.8341 x 10" 14 
0 
0 


(E.6) 


(E.7) 


Control actuation dynamics 
system of the form 


are modeled as a continuous linear time-invariant 


Xacf(0 — A ac (X acf(0 d* ^acf^aet(0 (E.8) 

u p(0 = C ac tx ac i(t^ -f" D act U a ct(t) (E.9) 


where x act (t) contain the actuator model states and u act (t) the actuator command 
inputs. For both flight conditions, we use the following actuator dynamic model, 


( ii(i) N 
* 2(0 
M0 

( M0 \ 

\ m ) 


-2 2 
0 -10 
0 0 

lOO’ 
0 0 1 


0 * 


' x,(0 ) 


0 o' 

0 


MO 

+ 

10 0 

-16 


s, 13(0 ) 


0 16 

xi(0 ' 


’ 0 

0 ' 

( M(0 \ 

X 2 (0 

+ 

0 

0 

( MO j 

M0 y 


L 




M(0 

M0 


(E.10) 
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where 



&T c ( 0 

MO 

h(t) 

Se(t) 


total thrust command (lbs) 
elevator position command (deg) 
total thrust output (lbs) 
elevator control position (deg) 


Table E.l: TSRV B-737 Trim Parameters at Flight Condition FLTl 


Parameter 

Value 

Altitude (ft) 

1500 

Gross Weight (lbs) 

80000 

Calibrated Airspeed (kts) 

120 

Center of Gravity Position 

0.2 

Flap Position (deg) 

40 

Flight-Path Angle (deg) 

-3 

Gear Position 

down 


Table E.2: TSRV B-737 Trim Parameters at Flight Condition FLT2 


Parameter 

Value 

Altitude (ft) 

25000 

Gross Weight (lbs) 

80000 

Calibrated Airspeed (kts) 

450 

Center of Gravity Position 

0.2 

Flap Position (deg) 

0 

Flight-Path Angle (deg) 

0 

Gear Position 

up 




Appendix F 

COMMAND PROCEDURES FOR DESIGN AND 

ANALYSIS 


Command procedures to set up the design synthesis models and perform robustness 
analysis for the AFCS controller design are written for the control system analysis and 
design software M AT LAB 1 . Command m-files for MAT LAB are used to develop the 
appropriate plant and controller models for use with the computer program SANDY ' 
and to perform multivariable robustness analysis. SANDY command files listed here 
contain the actual design values used in the final optimized AFCS designs. 

FA MAT LAB Command Files 

FA A Formulation of the Synthesis and Controller Models 

A M AT LAB function command file was used to compute state-space representations 
of the plant and controller models to be used as input data files to the computer 
program SANDY. A listing of this command file called msystem is given below. 

function msystem(f ilename ,f ilenamel , f ilename2) 

'/.// 

*/.// system(f ilename, filename 1 , f ilename2) 

'/,// computes the state-space representations of 

'/,// the plant (B-737) and controller (TECS) and closed loop system 
*1/ / for optimization to step commands in the constrained parameter 
*/•// optimization program SANDY. 

•/•// 

7 .// Data Input: A737W - Windaxis aircraft model (fps , rad, sec) 

7 .// B737W - Windaxis aircraft model (fps,rad,sec f lbs,deg) 

7.// VTBO - True airspeed (fps) 

7 // GRAVITY - Gravitational acceleration (fpss) 

Vf ALT - Altitude (ft) 

% U f LU - Ugust factor at altitude (fps) 


1 A trademark of Math Works, Inc. 
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X// LW - Wgust factor at altitude (fps) 

'/.// kgw - Aircraft gross weight (lbs) 

X// kti - thrust integral gain 

*/.// ktp - thrust proportional gain 

'/•// kei - elevator integral gain 

’/.// kep - elevator proportional gain 

X// kv - Outer loop gain 

'/,// kh - Outer loop altitude error gain 

X// Output Systems: 

’/.// 

X/ / ***** Formulate the Plant State-Space Matrices ***** 

7./ / System model : xpdot = Ap*xp + Bp*u + Gp*w 
7*// yp = Cp*xp + Dp*u + Omegap*w 

7.//where xp = [V(fps) .Alpha(deg) ,Q(deg) ,Theta(deg) ,Xdthl(lbs) t Xdth2 , Xde3(deg) , 
7#// Xgustv(fps) ,Xgustwl(fps) , Xgustw2,H(ft)] 

7.// u = [Dthc(lb) ,Dec(deg) f Vc*(fps)] 

7.// w = [Etav(fps) ,Etaw(fps) ,Gammac(deg) , Vdotc(f t/sec~2) , Vc(fps) t Hc(f t)] 

7.// y = [Gamma (rad) ,Vdot(g) ,Q (rad/sec) , Theta (rad) , (V-Vc) (fps) , 

7.// (Gamma-Gammac) (deg) , (H-Hc) (ft) ,Gammac(deg) ,Vdotc(ft/sec“2), 

7.// V(fps) ,Alpha(deg) ,Xdthl(lbs) ,Xde(deg) .Xgustv(fps) , 

7.// Xgustwl (fps) ,H(ft) ,Nz(g)] 

7, //where: x = [xIE,xIL,xVc] 

7.// uc = [Gamma(rad) , Vdot(g) f Q(rad) ,Theta(rad) # (V-Vc) (fps) , 

7./ / (H-Hc) (f t) , Gammac(deg) , Vdotc(ft/sec~2)] 

7// yc = [Dthc(lbs/lbs) ,Dec(rad) ,Vc*(fps)] 

7.Create the Plant Model and Controller Design Gains 
evaKC’load 9 , filename]); 

7.Conversion Constants 
r2d=180/pi; 7. //Radians to degrees 
d2r=pi/180; X //Degrees to radians 
XDryden Turbulence Filters 
alphau=vtbO/lul ; 

numv=[sqrt(2*alphau)] ;7. //numerator coefs. for u-gust filter 
denv=[l alphau] ; X //denominator coefs. for u-gust filter 

alphaw=vtbO/lw ; 

numw=sqrt(3*alphaw)* [1 alphaw/sqrt (3)] ;X //num. coefs. for w-gust filter 
denw=[l 2*.95*alphaw alphaw~2] ; X //den. coefs. for w-gust filter 

w2alpha=r2d/vtb0; X //conversion factor from w-gust to alpha-gust (deg) 

[av f bv ,cv ,dv] =tf 2ss (numv ,denv) ; 

[aw,bw,cw,dw] =tf 2s s (numw ,denw) ; 

[agust ,bgust ,cgust ,dgust]=append(av ,bv, cv,dv,aw,bw,cw,dw) ; 

X //State-Space Realization of u,w-gust Filters 
XGain Schedules 



149 


kgw=wg; '/,// Gross weight 

kcas = (200/vtb0) “2;'/, //elevator gain schedule 
'/.Actuator Dynamic Models 


CM 

CM 

1 

l l 

II 

0; 

*/, //actuator state model 

0 -10 

0; 

7. // 

x = [Xdthi ,Xdth2,Xdel] 

0 0 

-16] 

;7. // 

u = [Dthc (lbs/lbs) ,Dec(rad)] 

o 

i — i 

n 

0; 

V. 

// y = [Dth(lbs) ,De(deg)] 

10*kgw 

0; 




0 16*kcas*r2d] ; 
cact = [1 0 0; 

0 0 1 ]; 


dact = 0*eye(2); 

■/.Airplane Longitudinal Dynamics 
a737 = a737w ; 

a737 ( [2 3 4],:) = a737([2 3 4],:)*r2d; 

a737(:,[2 3 4]) = a737(:,[2 3 4])*d2r;7. //Convert A737w to fps and deg. 
b737 = b737w; 

b737([2 3 4],:) = b737([2 3 4],:)*r2d;7. //Convert B737w to fps and deg. 
a P - Ca737 b737*cact a737(:,[l 2])*[l,0;0,w2alpha]*cgust; 

zeros(3,4) aact zeros(3,3); 

zeros(3,4) zeros(3,3) agust] ; 

7. Addition of the Altitude State 
aplant = [ap zeros ( 10, i); 

0 -vtb0*d2r 0 vtb0*d2r 0000000]; 
bp = [zeros(4 , 3) ; 

[bact,zeros(3, 1)] ; 
zeros(4 ,3)] ; 

gp = [a737 ( : , [1 , 2] ) *dgust ; 
zeros(3 ,2) ; 


bplant 

•/.// 

7 .// 

'/.// 

'/.// 

cplant 


bgust 

0 0 ]; 

[bp gp zeros(ll ,4)] ; 

y = [Gamma (rad) , Vdot(g) , Q (rad/sec) .Theta (rad) , (V-Vc) (fps) , 

(Gamma-Gammac) (deg) , (H-Hc)(ft) .Gammac(deg) ,Vdotc(f t/sec*2) , 
V(fps) .Alpha(deg) ,Xdthl(lbs) ,Xde(deg) .Xgustv(fps) , 
Xgustwl(fps) ,H(ft) ,Kz(g)] 

[0 -d2r 0 d2r 0000000 
aplant (1 , :) /gravity 
0 0 d2r 00000000 
000 d2r 0000000 
10000000000 
0-1010000000 
00000000001 
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00 0 00000000 
00000000000 
10000000000 
01000000000 
00001000000 
00000010000 
[zeros (2 ,7) , cgust,zeros(2, 1)] 

00000000001 

vtb0/gravity*d2r*(aplant(4, :)-aplant(2, :)) ] ; 
dplant = [zeros(17,9)] ; 
dplant(5,3) = -1 ;dplant(E,8) = - 1 ; 
dplant(6,6)=-l ; 
dplant(7,9) = -1; 
dplant (8 ,6) = 1; 
dplant(9,7) = 1; 

dplant(17 , : ) = vtb0/gravity*(bplant(4, : )-bplant(2, : )) ; 

/•Save the Plant Synthesis Model into filenamel 

eval( [’save 1 , f il enamel , ’ aplant bplant cplant dplant']); 

’/.// 

I*/ f ***** Formulate the Controller State-Space Matrices ***** 

’/.// Initial 

’/,// System model xcdot = Ac*xc + Bc*uc 
*/.// yc = Cc*xc + Dc*uc 

’/.//where: x = [xIE.xIL.xVc] 

7.// uc = [Gamma(rad) ,Vdot(g) ,q(rad) .Theta(rad) , (V-Vc) (fps) , 

M t (H-Hc)(ft), Gammac(deg) , Vdotc(ft/sec“2)] 

’/.// yc = [Dthc(lbs/lbs) ,Dec(rad) ,Vc*(fps)] 

'/.// 

acont = [zeros(3,3)] ; 

bcont = [-1 -100 kv kh d2r 1/gravity 
1 -1 0 0 kv -kh -d2r 1/gravity 
00000 0 0 1 ]; 
ccont = [kti 0 0 

0 kei 0 
0 0 kl ]; 

dcont = [-ktp -ktp 00 0000 

kep -kep kq ktheta 000 0 

0 0 00 0 0 k2 0 ] ; 

’/.Save the Controller State model in filename2 
eval([’save ' ,filename2, * acont bcont ccont dcont']); 
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F.1.2 Multivariable Robustness Analysis 

Multivariable robustness measures described in Appendix C.3 involve computation 
of the minimum singular values of the return difference and inverse return difference 
transfer function matrices at different frequencies. MAT LAB function files were 
devised for computing the minimum singular values of these transfer function matrices 
from the state-space representation of the broken-loop systems. Osborne’s method 
[Ref.22] is used to improve the minimum singular values through diagonal scaling of 
the loop gain transfer function matrices. The MAT LAB command files svdrtd and 
svdirtd are user-defined functions for computing the minimum singular values of the 
return-difference and inverse-return-difference transfer function matrices respectively. 
Their listings are shown below. 

function [OMEGA , SIGMA] =s vdrtd( A , B,C,D) 
clg; 

J=sqrt(-1) ; 

[HIO,NIO]=size(D) ; 

[NS,NS]=size(A) ; 

DEL=0. 00001; 

NPTS=100 ; 

OMEGA =logspace(-2,2,NPTS) ; 

R = eye(HIO) ; 

RoldIMl=eye(NI0) ; 
for 1=1 :NPTS ; 

H = eye(NIO) + (C*inv(J*OMEGA(I)*eye(HS)-A)*B+D) ; 

MI = inv(M); 

Roldl = eye(NIO); 
rk = 2; 
k = 1; 

while abs(rk-l) > DEL; 
rsum = 0.0; 

for 1=1 :NI0;rsum=rsum+abs(MI(k,l))*2.O;end; 
rsum = rsum - abs(MI(k,k))“2.0; 
csum = 0.0; 

for 1=1 :lII0;csum=csum+abs(MI(l,k))*2.0;end; 
csum = csum - abs(MI(k,k))‘2.0; 
if rsum "= 0.0; 
if csum ”= 0.0; 

rk = (csum/rsum)*0.2S; 

R = eye(NIO) ; 

R(k,k) = rk; 



Roldl = RoldI*R; 
end; 

elseif rsum == 0.0; 
rk = 1.0; 

R = RoldIMl ; 

Roldl = R; 
elseif csum == 0.0; 
rk = 1.0; 

R = RoldIMl ; 

Roldl = R; 
end; 

H = R*M*inv(R); 

MI = inv(M) ; 
k=k+l ; 

if k > NI0; k = l;end; 
end; 

SIGMA(I)=min(svd(M) ) ; 

RoldIMl=RoldI ; 
end 

axis ( [-2, 2,0, 2] ) ; 
semilogx(0MEGA, SIGMA) ; 
xlabel( 'Frequency (Rad/Sec)’); 
ylabel( 'Sigma' ) ; 
grid; 

title( 'Minimum Singular Value of the Return-Difference Matrix’); 
rmin=min(SIGMA) 

’/, Multivariable stability margins 
MIM0GH=[-20*logl0(l+rmin) ,-20*logl0(l-rmin)] 

MIM0PM= [-acos ( 1-rmin* 2/2) * 180/pi , acos ( l-rmin*2/2) * 180/pi] 

function [OMEGA , SIGMA] =svdirtd( A , B,C,D) 
clg; 

J=sqrt(-1) ; 

[NI0 ,NI0] =size(D) ; 

[NS,KS]=size(A) ; 

DEL=0. 00001; 

HPTS=100; 

OMEGA =logspace(-2 ,2 ,NPTS) ; 

R = eye(HIO) ; 

RoldIMl=eye(NI0) ; 
for 1=1 :NPTS; 

H = eye(NIO)- inv(C*inv(J*OHEGA(I)*eye(KS)-A)*B+D) ; 

MI = inv(M) ; 

Roldl = eye(NIO) ; 
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rk = 2; 

* = l; 

while abs(rk-l) > DEL; 
rsum = 0.0; 

for 1=1 :NI0;rsum=rsum+abs(MI(k,l))‘2.0;end; 
rsum = rsum - abs(MI(k,k))*2.0; 
csum = 0.0; 

for 1=1 :NIO;csum=csum+abs(MI(l,k))“2.0;end; 
csum = csum - abs(MI(k,k))*2.0; 
if rsum "= 0.0; 
if csum "= 0.0; 

rk = (csum/rsum)*0.25; 

R = eye(NIO); 

R(k,k) = rk; 

Roldl = RoldI*R; 
end ; 

elseif rsum == 0.0; 
rk = 1.0; 

R = RoldIMl ; 

Roldl = R; 
elseif csum == 0.0; 
rk = 1.0; 

R = RoldIMl; 

Roldl = R; 
end; 

H = R*M*inv(R) ; 

MI = inv(M); 
k=k+l ; 

if k > NI0; k = l;end; 
end ; 

SIGMA(I)=min(svd(M) ) ; 

RoldIMl=RoldI ; 
end 

axis ([-2, 2, 0,2]) ; 
semilogx (OMEGA .SIGMA) ; 
xlabel( ’Frequency (Rad/Sec)’); 
y label ( ’Sigma’ ) ; 
grid; 

title( ’Minimum Singular: Value of the Inverse-Return-Difference Matrix’); 
rmin=min(SIGMA) 

‘/.Multi variable stability margins 
MIM0GM=[20*logl0(l-rmin) ,20*logl0(l+rmin)] 

MIM 0 PM= C-acos(l-rmin‘2/2)* 180/pi, acos(l-rrain* 2 / 2 ) ♦180/pi] 
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F.2 SANDY Design Data Files 

Data files for use with the computer program SANDY contain design information 
such as weighting matrices in the performance indices, linear constraints and nonlinear 
constraints. The SANDY data file for the optimized TECS inner-loop design in 
chapter 5 at flight condition FLTl is shown in Appendix G, while the one for flight 
condition FLT2 is shown in Appendix H. 



Appendix G 


DESIGN DATA FILE FOR TECS INNER-LOOP AT 
FLIGHT CONDITION FLT1 


Hf cmax' 

1 

1 

2000 

Tf f 

1 

1 

100 

Tfctor J 

1 

1 

2 

Npm’ 

1 

1 

4 

Wp’ 

4 

1 


1 

1 

, F > 40 10 

-0.04697 
-0.560599755740579 
8 . 5931e-14 
0.00423385423749794 
0.9995 
0 

-0.127046218711467 

-1.235 

0.000360974870088321 

-0.0125855088626523 

0 

0 

0 

0 

0 

2 
0 
0 
0 
0 
0 
0 


1 


0.298276769165831 

0.00040177 

-0.0234645506584104 

-0.0873187679779375 

0.008295 

-0.0441183231828685 

-0.0101291105986159 

-0.502 

0 

-0.219873991183662 

0 

0 

0 

0 

0 

0 

0 

0 

-10 

0 

0 

0 


1 


-4 . 00727596257898 e- 13 

0 

0.0531038896468495 

-0.7136 

-5.1027048276556e-06 
-0.0436213679933638 
-0.0251929542519023 
-0.0003233 
-1 . 13084680024971 
-0.0175300610836472 
1 
0 
0 
0 

-2 

0 

0 

0 

0 

0 

0 

0 
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-16 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-4 . 00727596257898e-13 

0 

0.0531038896468495 

-0.7136 

-5.1027048276556e-06 

-0.0436213679933638 

-0.0251929542519023 

-0.0003233 

-1.13084680024971 

-0.0175300610836472 

1 

0 

0 

0 

-2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.124782409638554 

0 

0 

0 

-0.0190695477624178 

0 


0 

0 

0 

0 

0 

0 

0 

-0.262375813333333 

0 

0 

0 

-0.04697 

-0.560599755740579 

8.5931e-14 

0.00423385423749794 

0.9995 

0 

-0.127046218711467 

-1.235 

0.000360974870088321 

-0.0125855088626523 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

-16 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

-0.124782409638554 

0 

0 

0 

-0.0190695477624178 

0 

0 

1 

0.298276769165831 

0.00040177 

-0.0234645506584104 

-0.0873187679779375 

0.008295 

-0.0441183231828685 

-0,0101291105986159 

-0.502 

0 

-0.219873991183662 

0 

0 

0 

0 

0 

0 

0 

0 

-10 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.262375813333333 

0 

0 
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0 

1 

0.298276769165831 

0.00040177 

-0.0234645506584104 

-0.0873187679779375 

0.008295 

-0.0441183231828685 

-0.0101291105986159 

-0.502 

0 

-0.219873991183662 

0 

0 

0 

0 

0 

0 

0 

0 

-10 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.262375813333333 

0 

0 

0 

-0.04697 

-0.560599755740579 
8 . 5931e-14 
0.00423385423749794 
0.9995 
0 

-0.127046218711467 

-1.235 


0 

0 

-4 . 00727596257898e-13 

0 

0.0531038896468495 
-0.7136 
-5. 1027048276556e-06 
-0.0436213679933638 
-0.0251929542519023 
-0.0003233 
-1.13084680024971 
-0.0175300610836472 
1 
0 
0 
0 

-2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.124782409638554 

0 

0 

0 

-0.0190695477624178 

0 

0 

1 

0.298276769165831 

0.00040177 

-0.0234645506584104 

-0.0873187679779375 

0.008295 

-0.0441183231828685 

-0.0101291105986159 

-0.502 


0 

-0.04697 

-0.560599755740579 
8 . 5931e-14 
0.00423385423749794 
0.9995 
0 

-0.127046218711467 

-1.235 

0.000360974870088321 

-0.0125855088626523 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

-16 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-4 . 00727596257898e-13 

0 

0.0531038896468495 
-0.7136 
-5 . 1027048276556e-06 
-0.0436213679933638 
-0.0251929542519023 
-0.0003233 
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0.000360974870088321 
-0 . 0125855088626S23 
0 
0 
0 
0 
0 
2 
0 
0 
0 
0 
0 
0 

-16 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


'G' 40 3 

0 

0 

0 

0 

0 

800000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

800000 


0 

-0.219873991183662 

0 

0 

0 

0 

0 

0 

0 

0 

-10 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.262375813333333 

0 

0 

0 


0 

0 

0 

0 

0 

0 

854.633071305922 

0 

0 

0 

0 

0 

0 

0 

0 

0 


-1.13084680024971 

-0.0175300610836472 

1 

0 

0 

0 

-2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.124782409638554 

0 

0 

0 

-0.0190695477624178 

0 

0 

1 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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'Gam 1 


0 

0 

0 

0 

0 

0 

0 

0 

0 

800000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

800000 

0 

0 

0 

0 

40 6 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

800000 

0 

0 

0 

0 

0 

0 

1 


854.633071305922 

0 

0 

0 

0 

0 

0 

0 

0 

0 

854.633071305922 

0 

0 

0 

0 

0 

0 

0 

0 

0 

854.633071305922 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

854.633071305922 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 



160 


o 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

800000 

0 

0 

0 

0 

0 

0 

1 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

800000 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

854.633071305922 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

854.633071305922 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 
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' Hs ' 


0 

0 

0 

0 

0 

0 

0 

0 

0 

800000 

0 

0 

0 

0 

0 

0 

1 

0 

0 

32 10 

0 

0.0174532925199433 

0 

0 

1 . 24526909962057e-14 
0 

0.00165021409716748 

0 

0 

0 

0 

0.0174532925199433 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

854.633071305922 

0 

0 

0 

0 

0 

0 

0 

-0.0174532925199433 

0 

0 

-0.00145960223741454 
-0.0174207506445177 
2 . 670323 18210068e- 15 
0.00013156787562144 
0.0174532925199433 
0 
0 
0 
0 
0 
1 
0 
0 
0 
0 
0 
0 
0 
0 
0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0.00926901085039872 

1.24850839030454e-05 

-0.000729165651286837 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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0 
0 
0 
0 

-0.0174532925199433 

0 

0 

-0.00145960223741454 
-0.0174207506445177 
2. 670323 182 10068 g- 15 
0.00013156787562144 
0.0174532925199433 
0 
0 
0 
0 
0 

1 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

0.00926901085039872 

1.24850839030454e-05 

-0.000729165651286837 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0.00926901085039872 

1.24850839030454e-05 

-0.000729165651286837 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0174532925199433 

0 

0 

-1 . 24526909962057e-14 

0 

0.00165021409716748 

0 

0 

0 

0 

0.0174532925199433 

0 


0 

0 

0 

0 

0.0174532925199433 

0 

0 

-1 .24526909962057e-14 

0 

0.00165021409716748 

0 

0 

0 

0 

0.0174532925199433 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.0174532925199433 

0 

0 

-0.00145960223741454 

-0.0174207506445177 

2.67032318210068e-15 

0.00013156787562144 

0.0174532925199433 

0 

0 

0 

0 

0 
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'Dsii 1 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 .0174532925199433 
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0 
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0 
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0 

0 
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0 
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0 

0 

0 

0 
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0 
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0 
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0 
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0 

0 

0 
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0 

-1 
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' Wspec * 


8 6 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

24 6 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

-i 

0 

0 

0 

-1 

0 

0 

0 

-1 

0 

0 

0 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

1 

0 

1 

0 

1 

1.2 

0 

1.2 

1.2 

0 

0.2 

0 

1 

0 

1 

0 

1 

0 
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’Alpha' 


’Sigmamax ’ 


’Zetamin’ 


’Q* 


o 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

4 1 

0 

0 

4 1 

0 

0 

4 1 

0.7 0 

0 

8 2 

1 
2 
0 
0 
0 
0 


0 

0 

0 

0 

0 

0.16 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 


0 

0 

0 

0 

0 


1 

0 

1 

0 

1 

0 

2 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

6.6 

1 

6.3 

1 

0 

1 

0 

1 

0 


0 


0 


0 

0 

0 

0 

0 
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, R* 12 3 

0 

0 

0 

10 

0 

0 

0 

0 

0 

0 

0 

0 

J nDircu* 4 1 

0 

2 

>Dircu> 2 3 

1 

2 

'nDircy 1 4 1 

0 

1 

’Dircy* 1 3 

1 0 5 . 1223e-l 

'A* 3 3 

0 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 
O.OOOOOOOOOOOOOOOD+OO 
’B 1 3 8 

-1 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 
1.745329251994330D-02 
1 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 
-1 .745329251994330D-02 
0 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 
'C' 3 3 

3.498961631412376D-01 
0 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 
’D’ 3 8 


0 

0 

0 

0 

0 

0 

0 

3 

0 

0 

0 

0 

0 


0 

0 

0 


0 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 
O.OOOOOOOOOOOOOOOD+OO 

-1 .OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 
3. 107520198881290D-02 
-1 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 
3. 107E20198881290D-02 
0 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 
1 . OOOOOOOOOOOOOOOD+OO 

0 . OOOOOOOOOOOOOOOD+OO 
2 . 403685052903295D+00 
0 . OOOOOOOOOOOOOOOD+OO 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0.000756 
6 . 735e-4 

0 


0 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 
O.OOOOOOOOOOOOOOOD+OO 

0 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 

0 . OOOOOOOOOOOOOOOD+OO 
-0 . OOOOOOOOOOOOOOOD+OO 

0 . OOOOOOOOOOOOOOOD+OO 
O.OOOOOOOOOOOOOOOD+OO 


0 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 
9 . 1941353801 13 152D-01 
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-9 . 6942329825 19398D-02 
0 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 
3 . 3759540980203 12D+ 00 
3. 124000504141355D+00 
0 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 
0. OOOOOOOOOOOOOOOD+OO 
-1 . 6332434490399 5 9D-01 


'nAid* 

’Aid* 

’nBid* 

'Bid* 

’nCid* 

'Cid* 

f nDid* 

'Did* 


3 7 

-1 1 

, Nclin > 

'nLincoef 1 


1 1 

0 

0 0 

i 1 

0 

0 0 

1 1 

3 

3 5 

1 

-5 

2 

5 

1 

1 1 

7 

7 5 

1 

-5 

2 

5 

1 

2 
-5 

3 

10 

1 


1 1 

2 

2 1 

2 

4 4 


-9 . 6942329825 19398D-02 
0 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 
-3 . 375954098020312D+00 
0 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 


1 

5 

1 

3 

0.85 


1 

5 

1 

2 

-5 

2 

5 

1 

2 

-10 


2 


0 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 

3 . 007108016451 16 1D+00 
0 . OOOOOOOOOOOOOOOD+OO 

0 . OOOOOOOOOOOOOOOD+OO 
0 . OOOOOOOOOOOOOOOD+OO 


1 

2 

-5 

3 

1 


1 

1 

-5 

1 

5 

1 

2 

-10 

4 

10 


Lincoef 1 
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1 

1 

1 

100 

1 

2 

’Linbnds 1 2 2 

0 

0 

* Nudnlc 1 1 1 

0 


100 

-1 

2 

2 

100 


0 


1 

100 

1 

1 

2 

0 



Appendix H 


DESIGN DATA FILE FOR TECS INNER-LOOP AT 
FLIGHT CONDITION FLT2 


* Nfcmar 1 

1 

1 

1000 

, Tf , 

1 

1 

100 

'Tfctor* 

1 

1 

2 

1 Npm * 

i 

1 


4 


4 1 

1 

1 

1 

1 

40 10 

-0.007055 

0.371755130674792 

-5.74387856831334e-13 

-0.561472420366576 

0.0004021 

0 

1.8924e-13 

-0.00655643326711755 

0.0320809006871565 

0.00799829068925111 

-0.007305211887918 

-1.032 

0.9976 

8.341e-15 

-5.835575143407440-07 

0 

-0.0605158023217175 

-0.00678896303972901 

-0.0890572497252435 

-0.0222034218500884 

0.00202368693240207 

-6.948 

-1.009 

0 

0.000360860278529295 

0 

-5.67990251047039 

0.00188067587892728 

-0.599583111522279 

-0.149485828502339 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-2 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-10 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



173 


-16 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

6 . 7438785683 1334e- 13 

0 

0.0320809006871566 

-1.032 

5 . 835575 14340744e-07 
-0.00678896303972901 
0.00202368693240207 
0 

-5.67990251047039 

-0.149485828502339 

1 

0 

0 

0 

-2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.431828457142857 

0 

0 

0 

-0.18647581639838 

0 


0 

0 

0 

0 

0 

0 

0 

-0.820474068571428 

0 

0 

0 

-0.007055 
-0.561472420366576 
1 . 8924e-13 
0.00799829068925111 
0.9976 
0 

-0.0890572497252435 

-6.948 

0 . 000360860278529295 
0.00188067587892728 
0 
0 
0 
0 
0 
2 
0 
0 
0 
0 
0 
0 

-16 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

-0.431828457142857 

0 

0 

0 

-0.18647581639838 

0 

0 

1 

0.371755130674792 

0.0004021 

-0.00655643326711755 

-0.007305211887918 
8 . 341e-15 
-0.0605158023217175 
-0.0222034218500884 
-1.009 
0 

-0.599583111522279 

0 

0 

0 

0 

0 

0 

0 

0 

-10 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.820474068571428 

0 

0 
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0 

1 

0.371755130674792 
0.0004021 
0.00655643326711755 
-0.007305211887918 
8 . 341e-15 
-0.0605158023217175 
-0.0222034218500884 
-1.009 
0 

-0.599583111522279 

0 

0 

0 

0 

0 

0 

0 

0 

-10 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.820474068571428 

0 

0 

0 

-0.007055 
-0.561472420366576 
1 . 8924e-13 
0.00799829068925111 
0.9976 
0 

-0.0890572497252435 

-6.948 


0 

0 

-5.743878568313346-13 

0 

0.0320809006871565 

-1.032 

-5 . 83557514340744e-07 
-0.00678896303972901 
0.00202368693240207 
0 

-5.67990251047039 

-0.149485828502339 

1 

0 

0 

0 

-2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.431828457142857 

0 

0 

0 

-0.18647581639838 

0 

0 

1 

0.371755130674792 

0.0004021 

-0.00655643326711755 

-0.007305211887918 

8.341e-15 

-0.0605158023217175 

-0.0222034218500884 

-1.009 


0 

-0.007055 

-0.561472420366576 

1.8924e-13 

0.00799829068925111 

0.9976 

0 

-0.0890572497252435 

-6.948 

0.000360860278529295 

0.00188067587892728 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

-16 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-5.743878568313346-13 

0 

0.0320809006871565 

-1.032 

-5.83557514340744e-07 

-0.00678896303972901 

0.00202368693240207 

0 
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0.000360860278529295 

0.00188067587892728 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

-16 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


'G 1 40 3 

0 
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0 

0 

0 

800000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

800000 


0 

-0.599583111522279 

0 

0 

0 

0 

0 

0 

0 

0 

-10 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.820474068571428 

0 

0 

0 


-5.67990251047039 

-0.149485828502339 

1 

0 

0 

0 

-2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.431828457142857 

0 

0 

0 

-0.18647581639838 

0 

0 

1 


0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

64.2101954444935 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 
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’Gam* 


0 

0 

0 

0 

0 

0 
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0 

0 

800000 
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0 

0 

0 

0 

0 

0 

0 

800000 

0 

0 

0 

0 

40 6 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

800000 

0 

0 

0 

0 

0 

0 

1 


64.2101954444935 

0 

0 

0 

0 

0 

0 

0 

0 

0 

64.2101954444935 

0 

0 

0 

0 

0 

0 

0 

0 

0 

64.2101954444935 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

64.2101954444935 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 



177 


0 0 
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